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ICKERS-ARMSTRONGS 


Head Office: VICKERS HOUSE BROADWAY LONDON: S-W:1 


‘—a section with undercuts 
in one operation?" 


DEFINITELY, sir, and a far more 

complicated section than this if 
you require it! The answer is 
“extrusion ”—any section which the 
designer’s ingenuity can devise can be 
produced in one operation, to final 
dimensions. 


Extrusions, forgings, castings— 
HIDUMINIUM Aluminium Alloys 
lend themselves readily to every 
manufacturing process. Their versatility 
makes them indispensable materials for post- 
war development in many branches of industry. 
If you use—or intend to use—aluminium 
alloys, the Development Department will be 
pleased to send you further details on request. 
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MERLIN 


There can be no doubt that the Merlin 61 is 
one of the outstanding achievements of the 
war. Rolls-Royce can be justly proud of it 


and the United Nations devoutly thankful. 
“ Aeroplane,’’ Nov. 20th 


ROLLS ROYCE 


AERO ENGINE 
For Speed and Reliability 


In the accelerated development in design 
and structure of scientific equipment and 
apparatus, when the elimination of dead- 
weight is a first essential, the new needs 
are adequately met by the employment of 
specialised Aluminium Alloys. In this de- 
velopment we shall keep pace, adapting 
our product to the requirements of its users, 
our system’ of production being standard- 
ised only for quality. For we produce 
Aluminium Alloys for every conceivable 


purpose. *Tomorrow it will be for instruments 


of peaceful discovery. 


SLOUGH: BUCK 
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AIRCRAFT LIMITED 
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T is not mere facetiousness to ; 
drawa parallel between the work devic 
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Should you desire, they will gladly work with 
you, or train your own staff in using Linatex. at all 
WILKINSON RUBBER LINATEX LTD., which 
GREAT TOWER ST., E.C.3 A Basie Material for i0 the 
Telephone: Mansion House 4333 industr 
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PRECISION 
in design, manufacture 
and testing: the basis 
of Lockheed success. 


AUTOMOTIVE PRODUCTS 


COMPANY LIMITED 


HYDRAULIC 


Balanced Judgment 


HIS is the, damping cylinder 
of a modern prismatic direct 

reading precision balance. The 
device consists of a metal cylinder 
inside which is a horizontal plate ; this 
plate is attached to the beam of the 
balance and floats free of the internal 
wall. It acts with a piston motion, 
thereby applying an air brake to the 
oscillation of the beam and bringing 
it to rest almost instantaneously. 
This simple mechanism is a vital 
invention in the development of 
the modern prismatic balance—the 
speed and accuracy of which is 
essential to the Scientific Control 
of Magnuminium magnesium alloys. 
It is the application of this control 
at all critical stages of manufacture 
Which maintains MAGNUMINIUM 


of the light alloy M AGNUMIN | UM 
The Lightest Bructural. Metal 
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Owing to the large demand, we are at present greatly 
restricted as regards the purposes for which this steel! 
can be 
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Specialists in 
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VOKES AIR FILTER 
FOR AERO-ENGINES. 


VOKES 


he the gallant forces of another age, charging with machine- 
like precision, with a Death or Glory spirit, carving their way to 
Victory! Their traditions, their indomitable spirit lives on, inspiring 
the mighty mechanised units of the British Army to-day. 

To ensure that engines develop the maximum efficiency it is essential 
that they are supplied with air, oil and fuel oil free from impurities. 
This is the job of work undertaken by Vokes Filters—not partially, but 
completely.. Vokes Air, Oil and Fuel Oil Filters are used in factories, 
offices and public buildings; by the manufacturers of aero engines, 
Diesel engines, petrol engines and power plants; in the engines of cars, 
lorries and tractors. Vokes Filters are made by engineers for engineers, 
and are used in all forms of industry. You can put any filtration 
problem to Vokes experts with the definite assurance of a successful 
solution, 


Vokes Filters with their 99%, filtration efficiency are not only approved 
and adopted by British and Dominion Governments, but are in use in 
the world’s biggest engineering shops. There is a Vokes Filter for all 
filtration problems. Consult the firm who specialize in this work. 
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“Yos, son, 

in 20 years time 
you'll be thanking me 
for that roof” 


Confidence in the future is a national duty. 


Confidence in the future of Cellactite is 
merely a question of accepting the experience 
of others who installed this permanent pro- 
tection years ago. 
Cellactite Roof (with Cellactite ‘Medway ’ 
Ventilators if needed) will resist wind, rain, 
frost, heat and even corrosive fumes as long 
as the building it covers. Under control 
direction the distribution of CELLACTITE is 
at present for essential purposes only. May 
we quote or send you full particulars ? 


(1) Bitumen-enveloped, corru- . 
gated steel sheet on to which are 
bonded CELLACTITE Asbestos- 
Asphalt felts conferring (2) 
impermeability from without, 

and (3) special acid-proof and 
non-drip properties. 


CELLACTITE & BRITISH URALITE LTD. 
Terminal House, 52 Grosvenor Gardens, 
London, S.W.1 
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Cellactite, Sowest. London 
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ACRYLIC RESIN SHEET 


* Does not shrink or warp with age 
Or exposure to varying atmospheric 
conditions. 

* Has very low moisture absorption 
and is unaffected by sca-water. 


‘PERSPEX’ is the perfect trans- 
parent material for aircraft 


1.C.l (PLASTICS) LTD. 


(A subsidiary company 
Imperial Chemical Industries Ltd.) 


Sales Offices at: Mill Hill, London, N.W.7; Oldbury, 
near Birmingham; Alderley Edge, Cheshire ; Bristol ; 
York; Newcastle-on-Tyne ; Glasgow ; Belfast; Dublin. 


THE MODEL 7 46-Range Universal 
AVOMETER 


Electrical Measuring Instrument 


A self-contained, precision moving-coil instrument, 
conforming to B.S. Ist Grade accuracy require- 
ments. Has 46 ranges, poidieg for measuring 
A.C. & D.C. volts, A.C. &.DC. amperes, res'stance, 
capacity, audio-frequency power output and de- 
cibels. Direct readings. No external shunts or 
series resistances. Provided with automatic com- 
pensation for errors arising from variations in 
temperature, and is protected by an automatic 
cut-out against damage through overload. 
Some delay in delivery of Trade orders is inevitable, 
but we shall continue to do our best to fulfil your 
requirements as promptly as possible. 
Sole Proprietors and Manufacturers : 
AUTOMATIC COIL WINDER & ELECTRICAL 
EQUIPMENT CO., LTD 
Winder House, Douglas St, Londen, S.W. 4. 
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British Mercantile Air ae. 
Fast cargo planes delivering 
necessities to war stricken 
countries; helping reconstruction; 
establishing new markets for 
Britain. Planes of revolutionary 
design propelled by the Ace of 
Blades—Weybridge Blades ! 
Wooden propeller blades of proved 
efficiency and superiority that do 
not fatigue like metal. Blades that 
are cheaper and easier to repair. 


WEYBRIDGE 
BLADES 


THE ACE OF BLADES 
FOR ALL AIRCRAFT PROPELLERS 


THE AIRSCREW CO., LTD. 
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WINGS FOR VICTORY 
WEEK - 


A great week ...a grand week... a glorious 
week... is coming your way— soon! Between 
the beginning of March and the end of June, 
every town and village will be holding a splendid 
savings week, their Wings for Victory Week. 
This is your business! It’s an R.A.F. Week; and 
that means it’s your week. The vital war work 
you’re doing is backing up the men who fly — 
wonderfully. But back ’em up still more— by 
seeing that your town gets its Victory Wings. 
Local Committees all over the country are 
humming with activity — arranging their savings 
targets, their opening ceremonies and a hundred 
and one other details. You’ll help your town, 
won’t. you? You’ve probably very little time — 
but you may have time to help in some way or 
another. In a very real way this is your Week. 
So whatever you do— give it your active support! 
See that your Week is a real success both with 
the savers and with the savings. Roll out the 
savings !| The sky’s the limit ! 


Go ald out for 
VICTORY WINGS 


Issued by the National Savings Committe 
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| Branch of Hawker Siddeley Aircraft Ltd. 
The brothers Montgolfier are credited with 

the invention of the ‘hot air’’ balloon, and 

their many experiments led eventually to a | 
flight of one and a half miles, carried out in | 
1783. Later balloonists attempted to sustain | 
flight by ‘taking. the fire with them,” and 
though many disasters occurred their 
perseverance added much to man’s store of 
positive knowledge. A leading source of 
positive knowledge to the pilot of to-day 
is provided by Weston Aircraft Instruments, 
which so faithfully guide and inform him 
from take-off to landing. 
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MONTHLY NOTICES 
FEBRUARY, 1043. 


Annuai Ceneral Meeting. 

Notice is hereby given that the Annual General Meeting of the Royal Aero- 
nautical Society, with which is incorporated the Institution of Aeronautical 
Engineers, will be held on Saturday, March 27th, at 3.0 p.m. in the offices of 
the Society, 4, Hamilton Place, London, W.1. 


AGENDA, 

1. To read the Notice convening the Meeting. 

2. To receive and deliberate upon the Report of the Council on the state 
of the Society and the Balance Sheets of Aerial Science Limited and 
Aeronautical Trusts Limited for the year ended December 31st, 1942. 

3. To receive the nominations for Council for the years 1943-1945. 

4. To announce the list of Fellows elected by the Council in accordance 
with Rule 4. 


5. To elect the Auditors for the ensuing year for Aerial Science Limited 
and Aeronautical Trusts Limited. 


6. Any other business. 
By Order of the Council, 
J. LAURENCE PRITCHARD, 
Secretary. 
Light refreshments will be served after the Meeting. 


Council Election. 


Attention is drawn to the following Rules for election of Council. Nominations 
for Council should be received by March 6th, 1943. Nomination forms may be 
obtained from the Secretary. Attention is drawn to possible delays in the post 
at the present time. 

RULE 17.—One half of the Council (excluding the President and the immediate 
Past-President if a member of Council) shall retire annually. The 
members who shall retire shall be those longest in office, except 
as provided in Rule 77. Retiring members of Council who have 
served for two terms in succession (four years) shall not be 
eligible for re-election until the next annual election, when they 
will be eligible. 

RULE 72.—The composition of the Council shall be as follows :— 

(a) Not less than sixteen shall be elected from the technical or scientific 
grades of the R.Ae.S., that is to say, from among the Fellows, 
Associate Fellows, Associates and Graduates. 

(b) Of these sixteen, four at least shall be Fellows. 
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RULE 73.—-Nominations of candidates for election to the Council must be 
received by the Secretary not less than twenty-one days before 
the Annual General Meeting, with an intimation in writing by the 


Candidates that they are willing to serve. Nominations must be 


signed by one supporter and two seconders, whe must be each 
entitled to vote in the R.Ae.S. 


New Year Honours. 

The following members were included in the New Year Honours for January, 
The Council wish to congratulate them on their well deserved Honours :— 


Acting Air Marshal J. E. A. Baldwin, C.B., C.B.E., D.S.O., O.B.E., 
A.D.C., R.A.F. 


C.B.E. 
Group Capt. P. W. 
C.B.E. (Military Division). 


Bulman, M.C., A.F.C., F.R.Ae.S. 


Acting Air Cmdre. K. L. Boswell, C.B.E., \.F.R.Ae.S. 


O:B.E. 


W/Cmdr. E. A. Gibson, A.M.Inst.C.E., A.F.R.Ae.S. 
E. T. Jones, A.R.Ae.S. 

Captain A. C. P. Johnstone, A.R.Ae.S. 

M.B.E. 


F, C. Bray, A.R.Ae.S. 


Fit. Lt. J. Healey, A.R.Ae.S. 
D. A. G. Jones, A.F.R.Ae.S. 


Fit. Capt. A. 


Lingard, A.M.Inst.C.E., A.F.R.Ae.S. 


W. A. Needs, A.R.Ae.S. 
Squadron Leader R. Alexander, A.F.R.Ae.S. 


Associate Fellowship Examination, December, 1942. 


kK. 


Beard 


IF. Brading 


Herbert-Caesari 


IF. Carr 

J. G. Carson 
B. Carter 

J. Chalwin 
Ik. Cowan ... 
A. Curry 

J. Detmold 
QO. E. Eden 


Theory of Internal Combustion 
Place) . 

Applied Mathematics. 

Design (Aero Engines). 

Theory of Machines (First Place). 

Applied Mathematics. 

Aerodynamics. 

Design (Aircraft). 

Applied Mathematics. 

Pure Mathematics. 

Applied Mathematics. 

Aerodynamics. 

Aircraft Materials. 

Design (.\ircraft). 

Pure Mathematics. 

Design (ero Engines). 

Aircraft Materials (First Place). 

Design (ero Engines). 

Theory of Machines. 

Applied Mathematics. 


Engines 


The following candidates were successful in the December, 1942, Associate 
Fellowship examinations :— 


(First 
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M. Ellins ... Design (Aircraft). 
Pare ... Design (Aero Engines). 


Theory of Machines. 
Applied Mathematics. 
G. L. Gunstone ... Design (Aero Engines)—First Place. 
Theory of Machines. 
Applied Mathematics. 


1, A. King... ... Applied Mathematics. 
R. A. McGill ... ... Strength of Aeronautical Materials (First Place). 
G, R. Matheson .... Theory of Machines. 
F. Milligan ... Design (Aero Engines). 
Theory of Machines. 
Applied Mathematics. 
B. Moss ita ... Applied Mathematics. 
W. A. E. Puntis  ... Pure Mathematics. 
LL. ot.. Leger .... Pure Mathematics. 
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R. D. Williams ... Design (Aircraft) —First Place. 

V. Willis ve ... Theory of Internal Combustion Engines. 
Applied Mathematics (First Place). 

A. J. Willmore ... Aireraft Instruments—Theory and Design. 


Baden-Powell Memorial Prize. 

The Baden-Powell Memorial Prize has been awarded to Mr. G. L. Gunstone, 
who was the best Student in the December, 1942, Associate Fellowship Examina- 
tions. 


Associate Fellowship Examination, May, 1943. 
Candidates for the May, 1943, Associate Fellowship Examination are reminded 
that entries must be in by March 31st, 1943. 


Craduates’ and Students’ Section. 

Members of the Graduates’ and Students’ Section are invited to the following 
lectures arranged by the Graduates’ and Students’ Section of the Institution of 
Automobile Engineers, 12, Hobart Place, S.W.1:— 

Sunday, February 14th, 1943, at 2.30 p.m. Paper by Mr. E. S. Grillo on 
Drawing Office Practice.” 

Sunday, March 14th, 1943, at 3.0 p.m. Paper by Mr. J. H. Millward on 
‘Gearbox.’ 

Sunday, April 11th, 1943, at 3.0 p.m. Informal Talk by Mr. L. Mantell on 
‘* The Practical Aspects of Thermal Efficiency.” 


Sunday, May 16th, 1943, at 3.0 p.m. Informal Talk by Mr. D. H. Bramley 
on ‘* Works Organisation.”’ 
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Election of Members. 


The following members were recently elected :— 


Associate Fellows.—Eric Thomas Charles Frost (from Student), Alfred 
Henry Higgs (from Graduate), Herbert Frederick Leonard Jenkins, 


Louis William Rosenthal. 


Associates.—John Sheppard Cumming, Shirley Donald Dickins (from 


Graduate), John Kenneth McKenzie. 


Graduates.—Douglas Edward Ackland, Hanus Bloch (from Student), 
Norman Alexander John Harry, Jack Joseph Molins. 


Student.—Henry Lambert Garbutt Sunley. 


Journal Binding. 


The following is a list of revised prices for binding the Society's Journals :— 


1942 Volume = 7S 


1941 Volume, including Library Catalogue 8s. 
1941 Volume, excluding Library Catalogue 7s. 
Library catalogue separately bound... 
Previous Volumes of Journal ... 


od. 
6d. 
od. 
od. 


Journals for binding should be sent direct to the Lewes Press, Limited, Friars 


Walk, Lewes, Sussex, and the remittance to this office. 


Advertisements. 


The fact that goods made of raw materials in short supply owing to war 
conditions, are advertised in the Journal, should not be taken as an indication 


that they are necessarily available for export. 


Annual Subscriptions. 


Members are reminded that their subscriptions are due on January Ist, 1943. 


The following are the current rates :— 


Home. 

Founder Members ... 2 2 0 
Fellows 4 4 0 
Associate Fellows 3 3 0 
Associates 2 0 
Graduates (age 21-25) 2 2 0 
(age 26-28)... 212 6 
Companions .. 


* £1 1s. od. without JOURNAL. 


Abroad. 
£ d. 
2 2 6 
2 2 6 
2 2 © 
212 6 
© 


All those members who have relinquished their Journal have been notified of 


the reduction in their subscriptions due on January Ist, 1943. 


Additions to the Library. 


Pamphlets in italics with location reference following in brackets. 
Books marked * or ** may not be taken out on loan. 


*B.a.273.—Ministry of Aircraft Production Official Reports. 
Focke-Wulf F.W. 190 A-3. (13 photos attached). 


Enemy Aircraft 
Examination of 


B.M.W. 801/A.1 Engine (in Dornier 217/E.1 Twin-Engined Bomber). 
(48 Figs. and Photos attached.) (Y.3.ili, 11 and 12.) 
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B.f.85.—How Fast do Birds Fly? By Frank W. Lane. (Article in ‘‘ Country 
Life,’’ Nov. 7th, 1941.) (Y.10.b.7.) 

BB.b.102.—Flugzeugleichtmetallbau. Matthiesen and Co., Berlin. 1938. 
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September, 1942. Dominion Bureau of Statistics, Ottawa. 

D.b.182.— Wings for Transportation. (Recent Developments in Air Trans- 
portation Equipment.) By T. P. Wright. Smithsonian Institution, 
Washington, U.S.A. (Reprint No. 3,654.) 1942. (Y.30.ii.58.) 

E.b.93.—On the Elastic Stability of a Rectangular Plate when subjected to 
a Variable Edge Thrust. By D. M. A. Leggett. (Reprint from Pro- 
ceedings of the Cambridge Philosophical Society, Vol. 31, Aug., 1935.) 
Cambridge University Press. 1935. (¥.8.a.7a.) 

*EE.h.7 and 8.—Gas Turbines and Jet Propulsion for Aircraft. By G. Geoffrey 
Smith. Flight Publishing Co. 1943. 3/6. (Two copies.) 

G.e..A.85.—Guide to the Properties of Wrought Aluminium Alloys. Wrought 
Light Alloys Development Association. (Information Bulletin No. 2.) 
1942. 

*].e.24.—Table of Circular and Hyperbolic Sines and Cosines. Federal Works 
Agency, New York. 1939. 20/-. 

L.b.4o.—* Flight ’? Handbook. (3rd Ed.) Flight Publishing Co. 1943. 6/-. 

L.c.24.—Advanced Flying By Captain Norman Macmillan. Allen and 
Unwin. 1943. 2/6. 

L.d.84.—Air Navigation. British Empire Edition. (2nd Ed.) By P. H. V. 
Weems. McGraw Hill. 1942. 35/-. 

L.d.85.—Astronomical Navigation. By W. M. Smart. Longmans, Green. 
1942. 5/- (2nd copy.) 

L.e.10.—An Introduction to the Study of Map Projection. By James 
Mainwaring. Macmillan and Co. 1943. 5/-. 

M.a.25.—Elementary Electricity for Radio Students. By W. E. Flood. 
Longmans, Green. 1942. 1/-. (PM.1.a.5.) 

M.b.47.—The Morse Code: Learning and Practice. By R. G. Shackel. 
(2nd Ed.) Longmans, Green. 1943. 1/-. (PM.2.c.7.) 

N.a.76.—The New Frontiers in the Atom. By Ernest O. Lawrence. 
Smithsonian Institution, Washington, U.S.A. (Reprint No. 3,674.) 
1942. (Y.30.1i.56.) 

N.a.77.—Mathematics and the Sciences. By J. W. Lasley. Smithsonian 
Institution, Washington, U.S.A. (Reprint No. 3,656.) 1942. (Y¥.go.1i.57-) 

*P.c.38.—Abridgments of Patent Specifications: Class 4—Aeronautics. 
(1899-1904.) Patent Office. 

**R.b.235.—Letter from Mr. W. Le Feuvre, President of the Balloon Club, 

electing Mr. Percival Spencer to membership of the Club. 1888. 

of **R.£.82.—Aeronautical Medals and Tokens. By Dr. F. John Poynton. 
Typescript. 1942. 

R.f.83 and 84.—The 1914 Tests of the Langley “‘ Aerodrome.’ By C. G. 
Abbot. Smithsonian Institution, Washington, U.S.A. (Miscellaneous 
Collections, Vol. 103, No 8.) 1942. (Y.30.i1) (Two copies.) 

S.d.80.—Sea Power. By ‘‘ T.124.’’ Jonathan Cape. 1940. 


aft S.d.81.—Italy’s Prophet of Air War Havoc. (General Douhet.) Stencil of 
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STRESSES IN SPACE-CURVED RINGS REINFORCING THE EDGES 
OF CUT-OUTS IN. MONOCOQUE FUSELAGES. 


By N. J. Horr, Ph.D., A.F.R.Ae.S., A.F.I.Ae.S.* 


1. INTRODUCTION. 


The object of this paper is to contribute one step toward the knowledge of 
the stresses in the neighbourhood of openings in monocoque aircraft fuselages 
by developing a procedure for the calculation of the stresses in spaced-curved 
rings used for reinforcing the edges of openings. The general appearance of 
such an opening is shown in the perspective view of Fig. 1. The procedure 
presented here presupposes a knowledge of the forces which act upon the ring. 
These forces depend upon the relative distortions of the reinforcing ring and of 
the rest of the fuselage, and the exact determination of their magnitude under 
different conditions of loading of the fuselage presents great difficulties. 

Such a rigorous calculation of the forces is not attempted in this paper. It is 
rather suggested to make the assumption that the elastic properties of the ring 
are identical with those of the cut-out portion. On this assumption the loading 
of the ring consists of the stresses which would prevail in the fuselage at the 
location of the ring, had no opening been provided. In reality the distortions 
of the reinforcing ring are, as a rule, greater than those of the corresponding 
portion of the undisturbed fuselage, and consequently the forces acting upon the 
ring are usually smaller than those calculated on the assumption suggested. In 
general, therefore, the assumption is on the safe side from the standpoint of the 
necessary strength of the ring. 

Characteristic features of the procedure suggested are the use of two 
systems of co-ordinates, one fixed to the fuselage, the other travelling along the 
centre line of the space-curved bar; the determination of all geometric and 
mechanical properties of the bar as functions of the length of its centre line; 
the use of vectors for the representation of directions, displacements, rotations, 
forces, and moments; and the consideration of the effect upon the moments of a 
non-uniform warping of the cross section of the space-curved bar. 

Although the calculations according to the procedure suggested are basically 
simple, they are rather laborious. For this reason a detailed investigation is 
carried out in order to determine the effect of different orientations of the principal 
axes of inertia of the cross sections, of different values of the ratios of flexural 
and torsional rigidities, and of a non-uniform warping of the cross sections. 
From these investigations conciusions are drawn as to when it is permissible 
to use simplified procedures. 

After the major part of the calculations presented here was finished, the writer’s 
attention was called to a paper by K. Marguerre. (Ref. 3.) In it a complete 
analytical solution of the moment distribution in one particular space-curved ring 
is given for some important loading conditions. In cases when the geometric 
and mechanical properties of the space-curved ring resemble closely enough those 
of the ring investigated by Marguerre, it is of obvious advantage to make use of 
his results rather than to carry out the calculations suggested in the present paper. 
The magnitude of possible deviations from Marguerre’s results is, therefore, 
investigated in some detail. 


* The Polytechnic Institute of Brooklyn, New York. The major part of the paper is taken 
from a dissertation submitted to Stanford University, California, in partial fulfilment 
of the requirements for the degree of Ph.D. The author acknowledges his indebtedness 
to Professor S. Timoshenko for the criticism, advice, and suggestions received from him. 
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In the Appendix the theory of the torsion of thin-walled open sections is 
presented. The theory was first given by Wagner (Ref. 7) and again discussed 
by Kappus (Ref. 2) in papers written in German, and published in English 
translation by the N.A.C.A. Although these N.A.C.A. publications are easily 
available, their lack of clarity makes a new detailed derivation appear useful. 


2. GEOMETRIC AND MECHANICAL PROPERTIES OF THE SPACE-CURVED Bar. 

The shape of a space-curved bar reinforcing the edges of an opening in a 
monocoque fuselage is seldom, if ever, known analytically. It is usually given 
by two or three normal projections in a layout drawing. If one establishes a 
fixed system of rectangular co-ordinates x, y, z, corresponding values of these 
co-ordinates can easily be determined for a number of points of the centre line 
of the bar. From the differences of these co-ordinates for adjacent points the 
length of an element of the are of the centre line of the bar, 

As = (Az? + Ay? + Az?):, (1) 
can be determined with any required degree of accuracy by spacing the points 
sufficiently closely. With a reference point s=o established on the centre line 
of the bar, the distance (measured along the centre line) of any point of the 
centre line from the reference point can be obtained by adding up the elementary 
lengths As. Thus the connection between the parameter s and the co-ordinates 
L, y, 2, is known, and may be presented in tabular form or by three graphical 
plots of the functions x=2 (s), y=y(s), and z=z (s), respectively. 

In the following calculations use is made of the unit tangent vector T. This 
is defined as a vector of unit length having the direction of the tangent to the 
curve and pointing toward increasing values of the parameter s. Its three 
components ¢t,, t,, f, in the direction of the fixed axes of co-ordinates are given 
by the equations 

where the dash ' denotes a differentiation with respect to s. If the functions 
xz(s), y(s), and z(s) are not known analytically, the components of the unit 
tangent vector may be obtained by graphical differentiation. For instance, the 
value of t, at any point o is equal to the slope of the curve a=- (s) at the point 
az=o. ‘Thus the values of the three components of the unit tangent vector can 
readily be determined and tabulated. 

In addition, the components of two more unit vectors are determined and 
tabulated. These are the vectors P, and P,, defined as vectors of unit length 
having the directions of the two principal axes of inertia of the section. Their 
sense is so chosen that the three vectors T, P,, P, correspond to a right-hand 
system. The directions of P, and P, can be taken from the layout drawings. 
If the shape of the surface of the monocoque fuselage is known analytically, the 
directions of P, and P, can often be calculated upon the assumptions discussed 
Secs; 

In the procedure outlined in the present paper two systems of co-ordinates are 
used simultaneously. In the fixed system the direction of the y-axis is taken 
as the direction of the axis of the horizontally disposed fuselage. The axes 
az and z are in a plane perpendicular to y, z is horizontal, z is vertical, and the 
positive sense is determined so that the axes correspond to a right-hand system. 
The travelling system is defined by the three mutually perpendicular unit vectors 
T, P,, and P,. The angles between the co-ordinate axes of the travelling system 
and those of the fixed system are, in general, different at each point of the centre 
line of the space-curved bar. 

When the cross-sectional dimensions of the space-curved bar are large as 
compared with the radius of curvature of the centre line of the bar, the stress 
distribution due to bending is not linear over the cross section. In such a case 
two other vectors may be chosen in preference to P, and P, to establish the 
travelling system of vectors in conjunction with the tangent vector. These two 
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new vectors are the unit principal normal vector and the unit binormal vector 
of the centre line of the bar. The advantage of such a choice is due to the fact 
that any very small element of a space curve may be considered a plane curve, 
the plane of which passes through the tangent vector and the principal normal 
vector. This plane is called the osculating plane. If the moment acting in any 
section of the space-curved bar is resolved into the directions of the principal 
normal and binormal, the distortions of the element of the bar may be calculated 
with the aid of the theory of bending of plane curved bars as presented in text- 
books on strength of materials. Although one is seldom compelled to make use 
of such a refinement in the calculations, the pertinent formule of differential 
geometry are given here :— 
The components of the unit principal normal vector P are 
Px=px" (8) py=py" (8) (8) (3) 
where p is the radius of curvature of the centre line of the space-curved bar. 
(1/p) is given by the equation 
(1/p) = (a"? + + (4) 
Since the binormal vector B is a vector of unit length perpendicular to both 
the tangent vector and the principal normal vector, it can best be both represented 
and computed as the vector product of 7 and P, that is 


j | 
B=) t, (5) 
Px Py Pz | 


where i, j, k are unit vectors parallel to the co-ordinates 2, y, z._ In all the above 
formule each dash denotes one differentiation with respect to the parameter s. 

Obviously, all the above formule may be used for determining the unit vectors 
P and B graphically, if the centre line is not known as an analytic function of s. 
Unfortunately, however, the accuracy of two consecutive graphical differentiations 
is not too good. 

As mentioned in the Introduction, all distortions of the bar are represented 
by vectors. A rotation vector is defined as a vector the length of which repre- 
sents (to some convenient scale) the angle of rotation, the direction of which is 
parallel to the axis of rotation, and the sense of which, as shown by the arrow, 
corresponds to the direction of the advance of a right-hand screw which undergoes 
the same rotation. moment vector is defined in the usual way, its sense is also 
taken as corresponding to the right-hand screw rule. 


3. THE Space-CurVED RING oF SoLID SECTION. 

The space-curved ring corresponding to a reinforcement of the edges of the 
cut-out in a monocoque fuselage as shown in Fig. 1 is six times. statically 
indeterminate. If a complete cut is made through it at any point, the unknown 
stresses in the cut may be replaced by a direct force (in the direction of T), two 
shear forces (e.g., in the directions of P, and P,), and by three moments the 
vectors of which again may be taken parallel to T, P,, and P,, respectively. 
The first of these corresponds to the twisting moment, and the second and the 
third to the two bending moments about the two principal axes of inertia of the 
section, respectively. These three forces and three moments may be chosen for 
the six unknowns of the problem. 

In the simple structure of the cut ring it is possible to determine the forces 
and moments acting in every section under any known external loads. Similarly, 
the forces and moments in every section due to the six unknowns can be deter- 
mined. Six linear equations can be obtained then for the determination of the 
six unknowns if one makes use of the conditions of continuity. These require 
the vanishing of the six components of the resultant relative displacement of the 
two sides of the cut. The six components are the three relative displacements 
measured in the direction of the vectors 7, P, and P,, respectively, and the three 
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relative rotations about them. The relative displacements can be calculated trom 
integrals extended over the total length L of the centre line of the space-curved 
bar. The integrals can always be solved by numerical or graphical methods. 
In the case of curved bars the cross-sectional dimensions of which are small 
as compared to their length, the deformations caused by direct forces and by 
shear forces are, as a rule, much smaller than those caused by bending and 
twisting moments. It is common practice, therefore, to calculate only the latter 
and completely disregard the tormer. This practice will be followed in the 
present paper. Consideration of the distortions due to direct force and shear 
adds materially to the computational work, but otherwise presents no difficulties. 


/ \ 


Fic. 1. 


Perspective view of a cut-out in a monocoque fuselage. 


Fic: 2: 


Centre line of a space-curved bar reinforcing the edge of a cut-out. 


38 
j 
\ 
tl 
/ 
> 
= O tw 
x co 
Si 
\ \ 
\ M. \ 
» 6) 
Mi cur 
of 
"(S20 S=O | 
7 
For 
exte 


ow 


STRESSES IN SPACE-CURVED RINGS IN MONOCOQUE FUSELAGES. 
The equations for the calculation of the unknowns may be developed in the 
following way. In Fig. 2 let the internal moments in a section o be called 
positive, if their vectors point in the positive directions of the co-ordinates on 
the side of the section that belongs to the portion of the bar extending from 
s=o to s=oa. The unknown force components F,,, Pyo, F,,, and the unknown 
moment components Myo, Myo, Mz. at s=o are drawn in directions opposite to 
the positive x-, y-, 2- directions, since they are acting on the side of the cut 
opposite to that defined above. It is convenient to calculate the moments caused 
by the external forces acting upon the bar (not shown in the figure) with respect 
to the fixed system of co-ordinates 2, y, z rather than with respect to the travelling 
system ft. p,. p,. If these external moment components are identified by the 
subscript ¢. one may give the resultant moments acting in any section in the 
form :— 
M,=M,.+ M.. o% Foy 
My = Myo + Myo t+ Fact — (6) 
t+ Mio t+ — Fyot 
Before proceeding to the calculation of the distortions, one must resolve these 
moments in the directions of the travelling system of co-ordinates. It is known 
from vector calculus that the component of a vector in the direction of a unit 
vector is the scalar product of the two vectors. Hence the component of the 
moment .M (given in Eqs. 6 by its three components M,, M,, M,) in the direction 
of the tangent to the centre line of the bar is the scalar product of M and the 
unit tangent vector 7’ :— 
M.=M .T=M,t,+M,t,+Mga, . (7a) 
Similarly, the moments in the directions of the two principal axes of inertia 
of the section are :— 
M,=M .P,=M + ‘ ‘ 
The distortions of an element of the bar of a length ds can be given by the 
three angles of rotation dB,,, and df,., 
dB, =M,ds/GC 
dB,,=M,,ds (Kl, ‘ (8) 
where GC is the torsional rigidity of the has according to the Saint-Venant theory, 
I, and I, are the two principal moments of inertia of the section, and FE is Young’s 
modulus of the material. For the calculation of the relative distortions of the 
two sides of the section at s=o the rotation d@ (as given in Eqs. 8 by its three 
components df,, dB,,, and d&,,) must be resolved into the a-,y-, z- directions. 
Since the direction of the rotation vector df, is the same as that of the unit vector 
T, and since the 2-component of T is f,, the z-component of df, is the product 
t.dB,. Dealing with the other components in a similar way, one obtains 
dp, t.dB, + Poxd Boo 
dB, =tydBy + + Boo - : : (9) 
AB, = By + + Pox Boo 
The relative rotations of the two sides of the cut at s=o can now be calculated 
by integrating the elementary rotations over the closed centre line of the space- 
curved “bar. The integrals so obtained must vanish because of the continuity 
of the bar. Thus one ‘obtains the three equations 


odp, =O 
odB,=0 ; : (10) 
od8,=0 


The relative vertical displacement dw at s=o due to a rotation df is 


dw=adB, —ydB,. 


For reasons of continuity the integral of this expression must also vanish when 
extended over the closed centre line of the bar. From the requirements of a 
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vanishing of the relative displacements in the z- and y- directions, one obtains 
two more equations. Altogether one has 
oad B, 
oydB,—ozdB,=o0 }. ‘ (11) 

The integrals of the six equations Eqs. 10 and 11 can be calculated numerically, 
or graphically. All of them contain, in general, the six unknowns. Thus six 
linear equations are obtained which can always be solved for the six unknowns 
they contain. Substituting the values of the unknowns into Eqs. 7, one obtains 
the twisting and bending moments as functions of the parameter s, and thus 
the problem is solved. 

The procedure explained above, though always possible, is rather laborious. 
Fortunately, most space-curved bars used in aircraft are symmetrical with respect 
to one or more axes. When this is the case, the properties of symmetry may 
be used to decrease the number of the unknowns. 


PLAN View 


SECTION 


Fic. 3. 
Orthogonal projections of the centre line of the space-curved 
bar reinforcing the edge of the cut-out shown in Fig. 1. 


It should be noted that the above calculations are valid only if Hooke’s law 
is valid, and the distortions are small as compared with the cross-sectional 
dimensions. Moreover, the radius of curvature of the centre line of the bar must 
be large as compared with the cross-sectional dimensions. This last requirement, 
however, need not be fulfilled if the two principal axes of inertia of every cross 
section coincide with the principal normal, and the binormal of the curve, 
respectively, and the moments of inertia are replaced by their equivalents in the 
theory of plane-curved bars. It is further assumed that the shape of the cross 
section does not change under the loads (e.g., due to flattening or buckling). 


4. CALCULATION OF THE MoMENT DISTRIBUTION IN A DoUBLY SYMMETRICAL RING. 
As an example, the moment distribution in the doubly symmetrical space-curved 


bar, the centre line of which is shown in Fig. 3, is calculated. The curve may 
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STRESSES IN SPACE-CURVED RINGS IN MONOCOQUE FUSELAGES. 41 
be considered as representing the edges of a cut-out in the top portion of a 
smaller size monocoque fuselage. The calculations are carried out analytically 
rather than graphically in order to facilitate the investigation of the effects upon 
the moment distribution of changes in the values of the three rigidities GC, EI,,, 
and EI,,. 

The loading of the ring consists of two equal and opposite forces of a magnitude 
2P each acting parallel to the y-axis at the topmost points of the curve. In 
Fig. 4 one-half of the ring is shown together with the loading. The location 
of the origin of the fixed system of co-ordinates differs from that in Sec. 3. The 
effect of this change upon the equations of Sec. 3 is readily found. It is assumed 
arbitrarily that the directions of the principal axes of inertia of every section 
coincide with the principal normal, and with the binormal, respectively, of the 
curve. 


Fic. 4. 
Axonometric view of one-half of the space-curved bar of Fig. 3. 


The co-ordinates of the curve are given by the equations :— 
y= (a/3) [3 (42) 
a= —(a/3)72 [1-—(8/a)? 8? 


eS (s/a)S1 
where a is the length of one-quarter of the complete space-curved ring. The 
complete ring is obtained by mirroring the above-defined curve in planes parallel 
to the a- and z-, and to the y- and z-axes, respectively. 
The three components of the unit tangent vector T of the curve were calculated 
from Eqs. 2. They are 


The equations are valid if 


t.=(s/a)? 
v2 (8/a) [1—(8/a)?]} 
The principal normal vector P can be calculated from Eqs. 3. Its components 
are 
2 (8/a) [1 (8/a)?]? 
Py = — (8/a) [1 —(8/a)? : (14) 
Pz=1—2 (s/a)? 
The length of the radius of curvature is calculated with the aid of Eq. 4. It is 
p=(a! 2) [1—(s/a¥]} (15) 
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The binormal vector 4 is calculated from Eqs. 5. Its components are 
b,=1—(8/a)* 
by =(8/a)? (16) 
b,=— (s/a) [1 

Due to the double symmetry of both ring and loading, stresses and distortions 
are identical in the four quarters of the ring. It suffices, therefore, to consider 
only one quarter. If a cut is made both at s=o and s=da, in Fig. 4, the following 
force and moment components must vanish in the cuts because of the symmetry :— 

At s=o: the shear force parallel to 2, the shear force parallel to z, and the 
torque (moment vector parallel to y); 

At s=a: the shear force parallel to z; and the torque (moment vector parailel 
to x). 

Moreover, and still because of the symmetry, the tensile force at s=o (force 
parallel to y) is #, and the shear force parallel to y at s=a is also P. The values 
of some more force and moment components follow then from the requirements 
of static equilibrium : 

At s=o: the moment vector parallel to 2 is — Py 

At s=a: the tensile force (parallel to x) is zero, and the moment vector parallel 
to y IS zero. 

Thus the only two unknowns in the cuts are the moment vectors parallel to z. 
The two, however, are interdependent. If one denotes the unknown moment at 
s=o by M,,, the moment at s=a becomes M,,—Ps. With this notation one 
obtains for the moments at any point of the centre line of the ring :— 

M,=Pz 
M,=M,,-- Pz 
These equations correspond to Eqs. 6 of the preceding section. They contain a 
single unknown, namely \M,,, which must be calculated from the condition of 
continuity. Substituting the expressions of Eqs. 13-17 into Eqs. 7 and remem- 
bering that the vectors P, and P, are coincident with the vectors P and B, 
respectively, one obtains :— 
M,=—(Pa/3)¥ 2 [1 —(s/a)* Mov 2 (8/a) [1 —(s/a)?]}} | 
p= — (Pa/3) (8/a) [2-3 (s/a)?]+M,o — 2 (8/a)*} 
(Pa/3)¥/2 [1—2 (s/a)*] [1—(s/a)? 2 (8/a) 

The values of these twisting and bending moments were calculated in terms 
of Pa and M,,. They are plotted in Fig. 5, where the subscript ? refers to the 
former part of the moment, M to the latter. If one now calculates the relative 
distortions of the sections s=o and s=a, instead of extending the integration 
of the elementary distortions over the total length of the closed curve, the 
requirements set forth in Eqs. 10 and 11 must be revised. It is easy, however, 
to establish end conditions for the rotations, and to use them for the determination 
of the unknown moment M,,. The actual values of the rotations 48 must fulfill 
the following requirements of symmetry :— 

At s=o: AB,=o Af,=any value Af,=o 
At s=a: Af,=any value AB,=o Ap, =0 

Accordingly, the rotation about the z-axis must vanish at both ends of the 

quarter-curve. Consequently 


— 
ll 


(18) 


a 


Eq. 19 now replaces Eqs. 10 and 11. The definite integral on its left side can 
be calculated if one substitutes Eqs. 18 into Eqs. 8 (remembering that the 
subscripts p and ) replace the subscripts p, and p,, respectively), and the values 
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so obtained into the last of Eqs. 9. After integration between the given limits 
Eq. 19 becomes :— 


(1/G@C) [(— 1/18) Pa? + (4/15) M,.a] + (1/EI,) [(—1/12) Pa? + (7/15) M,.a] 


+(1/ETpy) [(1/18) Pa? + (4/15) M,a]=o . . (19a) 
The solution of Eq. 19a may be written in the form 
and the parameters g and h are defined as 
g=(GC/EI,) and h=(GC/EI,) .  (20b) 
O 
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Fig. 5. 
Moment distribution along the space-curved bar. 
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Calculated values of fy are plotted against g and h in Fig. 6. This figure 
permits an analysis of the effect of the torsional and flexural rigidities upon the 
redundant moment M,,._ This moment changes most rapidly with a varying E Ip, 
most slowly with a varying EJ,. A variable EI, corresponds to a vertical line in 
Fig. 6, a variable E I, toa horizontal line, and a variable torsional rigidity GC 
toa straight slanting line through the origin, if in each case the other two 
rigidities are held constant. Fig. 6 shows that the redundant moment is greater, 
the greater the flexural rigidity FJ, about the binormal, and the smaller the 
torsional rigidity GC. The flexural rigidity FJ, about the principal normal is of 
lesser importance. 


A 
10 
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Fic. 6 
Plot of fy of Eq. 20 against gy=CG/EI, and h=GC/EI,. 
The principal normal and the binormal are the principal 
axes of inertia. 


If the moment diagrams of Fig. 5 are also taken into consideration, it may be 
seen that the effect of the redundant moment M,, is to reduce the torsion of the 
bar. With a small value of GC and a high value of EJ, the resultant torque M, 
almost vanishes, but at the same time the maximum values of the bending 
moments are not increased. The shape of the moment diagrams, however, is 
greatly affected. On the other hand, the resultant moment diagrams are very 
ane affected by the redundant moment if the torsional rigidity and the ratio 

I,/I, are large. 

The details of the above conclusions naturally hold true only for the curve 
considered here, and for the assumed orientation of the principal axes of inertia. 
The effect of a change in the latter is investigated in the next section. It may 
be stated generally, however, that if the torsional rigidity GC is small, the 
redundant moment tends to cause a decrease in the values of the twisting moments 
so that the bar carries the external loads mainly by virtue of its bending rigidities. 
The same conclusion was drawn by Marguerre from his investigations. 


. Tue Most LikeLty DIRECTIONS OF THE PRINCIPAL AXES OF INERTIA. 
os many cases the section used for reinforcing the edges of a cut-out in a 
monocoque fuselage is intentionally given a permanent twist during the manu- 
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facturing process so that one of its plane surfaces is everywhere approximately 
tangential to the surface of the fuselage. In this way the section can easily be 
riveted to the covering of the fuselage. If the section is an J, or a channel, one 
principal axis of inertia is then coincident with the normal of the surface of the 
fuselage; if an angle, or a Z, the angle between the surface normal and the 
principal axis of inertia is constant. The direction of the surface normal, 
however, can be calculated if the shape of the surface is known analytically. 
It is shown in differential geometry that if the co-ordinates 2, y, z of the surface 
are known functions of two parameters u, v:— 
(u, v) 
y=y (1, 0) 
2=2 (u, v) 
and the differential coefficients of these functions exist and can be calculated with 
respect to both uw and v, and are represented by the symbols 
Ly=Or/OU Yy=Oy/OU = 02 / OU 
e,=02/d0  yy=dy/dv 2,=02/dv 
then the surface normal N is 


| } k | 

| | 
| &y Yy ay | 


This is not a unit vector. It can be transformed into one by dividing it by its 
length. 


If the surface of the fuselage is cylindrical, its parametric equation is easily 
found from the equation of the space-curved bar. If, for instance, the generatrices 
are parallel to the y-axis, and the centre line of the bar is given by the equations 

e=2(s) y=y(s) (8) 
then the parametric equation of the surface may be written in the form 
z=2(s) you (8) 
where s and v are the independent parameters. Thus the surface normal becomes 


] k 
N=| z! (s) fo) 2! (8) | 
I 
Accordingly the components of the surface normal are 
(8) 2, (2) 


The vector can be normalised by dividing it by its absolute magnitude. Hence, 
the components of the unit surface normal vector are 
n, = —2! (s)/[x!? (8) +2! (s)]} 
%, . (21a) 
n,=a! (8)/[a!? (s) +2! 

If the surface normal is one principal axis of inertia, the other is perpendicular 
to both the surface normal and the tangent of the curve. Consequently, the 
vector product of the unit tangent vector and the unit surface normal vector is a 
unit vector R which points in the direction of the second principal axis of inertia. 
It is given by the equation 


| j k 
R=) t, t, : (22) 
a. Ny n, 


6. THe EFFECT OF THE DIRECTIONS OF THE PRINCIPAL MOMENTS OF INERTIA UPON 
THE REDUNDANT MOMENT. 

With a view to investigating the effect of the directions of the principal moments 
of inertia upon the magnitude of the redundant moment M,,, the calculations of 
Sec. 4 are repeated here on the assumption that the directions of the principal 
axes of inertia are those of the vectors N and R defined in the preceding section. 
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In the case of the cylindrical surface defined by the curve given in Sec. 4 the 
components of the unit surface normal vector are :— 
Ny= — 
ny=O 
J 
Ihe components of the unit vector R are :— 
ry= —(8/a) [1 
=(s/a)[2 —(s/a)?}} (24) 
With these values the twisting and bending moments become :— 
M,= —(Pa/3) ¥ 2 (s/a)? 2 (8/a) [1 —(s/a)?]} 
M,=(Pa/3) [2—4 (s/a)? + (8/a)*] / a)? }} 
+ M,, (25) 
M,=(Pa/3)¥ 2 (s/a) [1 —(s/a)? 
— 2 [1 
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Fic. 7: 
Plot of fy of Eq. 26 against g=GC/EI, and h=GC/EI,. 
The surface normal is one principal axis of inertia, the other 
is perpendicular to the surface normal, and to the tangent 
to the centre line of the bar. 


Carrying out the calculations in the same way as in Sec. 4, one obtains 


where 
fy = [-0555 —-01896 g + .04673 h]/[.266+.244 9+-4854h] (26a) 
and the parameters g and h are the same as those defined in Eqs. 20). Calculated 
values of fy, are plotted against q and i in Fig. 7. 
Marguerre in his calculations (see Ref. 3) assumes that one of the principal 
axes of inertia of the section is always horizontal. The manufacturing of such 
a space-curved bar is comparatively easy, its attachment to the sheet covering 
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of the fuselage, however, is difficult. The moment distribution in the bar, the 
centre line of which was defined in Sec. 4, is now calculated on the assumption 
that one of its principal axes of inertia is horizontal. 

The vector pointing in the direction of the first principal axis of inertia must 
then be perpendicular to both a vertical unit vector and the tangent vector. 
Denoting this vector by H, one has 

] 


| (8/a) [1—(8/a)?}} | 


Consequently 
h,=(s/a)?-1 h,=(s/a)? h,=o 
The unit vector of the same direction has the components 
| he=[(8/a)? >—1]/[2 (s/a)*—2 
hy =(8/a)?/ [2 (s/a)*— 2 (s/ 
h,=o 


SO 


4 
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Fic. 8. 
Plot of fy of Eq. 30 against g=GC/EI, and h=GC/FEI,. 
One principal axis of inertia is horizontal, the other is 
perpendicular to the first principal axis, and to the tangent 
to the centre line of the bar. 


The unit vector N having the direction of the second principal moment of inertia 
of the section can be calculated as the vector product of the unit vectors H and T. 
Its components are 


(26) n,=¥v 2 (s/a)® [1 —(s/a)?]3/[2 (8/a)'—2 (s/a)? +1 

ny = 2 [(s/a)—(s8/a)®] [1 —(s/a)? ]?/[2 (s/a)*—2 (s/a)? +1 (28) 
26a) n,= —[2 (s/a)'—2 (s/a)? +1]! 
ated | The moment distribution is given by the equations 

M,= —(Pa/3)¥ 2 (s/a)? [1 —(s/a)? 2 (8/a) [1 —(8/a)?] 
cipal | [(s/a)?—1]°/?/[2 (s/a)* +1]3 (26) 
such | = — (Pa/3) (s/a)® (8/ay “T/T? (s/a)*—2 (s/a)?+1]} 


ring —M,, [2 (s/a)*— 2 (s/a)?+1]} 


(27) 
| 
i 
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For the unknown moment the expression is obtained 


where 
and the parameter g has the same meaning as before. It may be seen that in 
the present case the redundant moment is independent of h. Computed values 
of fy are plotted against g and hi in Fig. 8. 

It should be noted that the twisting moment M is given by the same expressions 
in all the three cases discussed, as may be seen from a comparison of Eqs. 18, 
25, and 29. The magnitude of fy, is, therefore, representative of the decrease 
in the values of the twisting moment due to the redundant moment M,,. A 
comparison of the Figs. 6-8 reveals that fy is independent of the directions ot 
the principal axes of inertia if g=h. On the other hand, the differences in the 
value of fy for different orientations of the principal moments of inertia are 
material, if one of the parameters g and hi is small, the other large. Consequently 
the values of the twisting moment depend considerably upon the orientation of 
the principal axes of inertia if the bending rigidity of the space-curved bar is 
materially different in the two principal directions, and if the torsional rigidity 
GC is not too small. 


7. NoN-UNIFORM WARPING. 

The calculation of the redundant moment M,, in Sec. 3 was carried out upon 
the assumption that the unit angle of twist of the curved bar was everywhere 
proportional to the torque, as expressed in the first of Eqs. 8. This is, however, 
true only if Saint-Venant’s assumptions concerning simple torsion, as stated in 
Sec. A.1 of the Appendix, are satisfied. If the torque varies along a bar of 
thin-walled open section, the angle of twist has to be calculated from the 
differential equation 

GC (dB,/ds) — EY (d*8,/ds*)=M, (31) 
This equation is derived in Sec. A.5. The symbol I’ denotes the warping constant 
which is defined in Eq. A.28. The meanings of the other symbols were previously 
explained. 

The seemingly inconsistent procedure of calculating the angle of twist of an 
element from a differential equation, while the angles of rotation due to bending 
are determined from the much simpler formule of Eqs. 8, has a simple physical 
reason. In bending, the well-established assumption that plane sections before 
bending remain plane after bending eliminates all interference between the elements 
of the bar during distortions. Thus the distortion of each element can be 
calculated individually. On the other hand, in torsion the effect of the warping 
of the originally plane cross sections is of primary importance and if it is not 
considered entirely erroneous results are obtained for both the angle of twist and 
the shearing stress, as was shown by Saint-Venant a century ago. Thus, if 
due to a variable torque, or to some restraint to warping, the warping is not 
uniform over the whole bar, the distortions of each element are influenced by 
the distortions of all the other elements. 

In the case of solid, or thick-walled closed sections the calculation of the effect 
of non-uniform warping is much more involved. Fortunately, with such sections 
the effect is negligibly small, except in the immediate proximity of concentrated 
loads and constraints. The angle of twist can be calculated, therefore, from 
the first of Eqs. 8 without undue error. For solid, or hollow circular sections 
the first of Eqs. 8 is rigorously correct, since such sections do not warp when 
twisted. 

The edge reinforcements most often encountered in monocoque aircraft, how- 
ever, are made of thin-walled open sections for which both terms on the left side 
of Eq. 31 are of equal importance. Even with sections of extremely thin wall 
thickness, for which EI’ may be much larger than GC, it is not advisable to 
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disregard the first term on the left-hand side of Eq. 31. If the term is disregarded, 
the equation is obviously incorrect in the case of a constant torque M,,. 
Eq. 31 may be given in the form 
where 
each dash ' denotes one differentiation with respect to s, and the symbol M, (s) 


is used in order to denote that the torque is a function of the parameter s. The 
solution of the reduced equation 


can be written as 
or 


The solution of the complete equation Eq. 31a is 


+ Bu! (34) 
where £,,' is a particular solution. This can be found by the operational method 
(K? — D?) B,,'=(1/ET) M, (s) 

Buy! =(1EV) [1 /(K?—D2)] M, (3) 
The latter equation may be transformed into 
Bi! = (1/2KET) { [1/(D+K)] —[1/(D—K)] } Me (8) 
From this form it is clear that 8,’ can be obtained as the sum of the solutions 
of two linear differential equations of the first order. These solutions are well 


known. Since only the particular solutions are sought here, the following 
expression is obtained :— 


= (1 /2KET) (s) ds — (s) ds (35) 


When this equation is substituted into Eq. 31a, an identity is obtained. Hence 
the differential equation Eq. 31a is solved for any type of distributed torque and 
the determination of the angle of twist is reduced to the solution of the two 
integrals in Eq. 35. ‘This can always be accomplished with the aid of numerical, 
or graphical methods. The two constants of integration in Eq. 33a (or in 
Eq. 330) must be determined from the boundary conditions. 

The redundant moment M,, can now be calculated according to the procedure 
discussed in Sec. 3, or in the case of double symmetry by that discussed in Sec. 4. 
The only difference is that in the calculations the values of £,! obtained from the 
differential equation Eq. 31 must, of course, be used in place of the values 
computed from the first of Eqs. 8. 


8. NUMERICAL EXAMPLE. 


In order to obtain a numerical appraisal of the effect of non-uniform warping 
upor M,,, the example of Sec. 4 is continued here. With the arbitrary assumption 
that 

K=1/a 
where a is again the (developed) length of one-quarter of the space-curved bar, 
Eq. 35 becomes 


Bp! =(a/2ET) (s) ds — em M, (s)ds] . 


After substituting the value of M,(s)=M, from the first of Eqs. 18 into 
Eq. 354, one has to integrate the terms containing Pa, and those containing 
M,. separately. The integration was carried out numerically. Its results for 
the ‘* moment ’’-terms (those containing M,,) are given in the table below :— 
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Writing the solution of the complete equation in the form 
=A cosh (s/a)+ B sinh (s/a) + (Bip’) (36) 
one can determine the values of the constants of integration trom the end condi- 
tions. Due to the double symmetry the warping must vanish at both end sections 
ef the bar. As may be seen from the derivations of the Appendix, this is the 
case when 8,/=o0. Hence the two end conditions give 
when s=o 
B,{=o0 when s=a 
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FIG. 9. 


Moment distribution 74=Miy/¥2M,, and unit angle of 
twist distribution (hq. 38). 
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Since (P,)')\, is zero when s=o, one obtains 
A=o 
For s=a the particular solution is 


(Bey!) = — -288 / 2M,,a?/2ET' 

Hence 
B=.288 / sinh (1)=.244 2M,,a?/2ET 

Substituting these values into Eq. 36, one may write 


(Bem 2M,,a7/2ET 


where 6, is a numerical factor. Its values are plotted against (s/a) in Fig. 3. 


(38) 


The figure also contains calculated values of 7, defined as the factor by which 
¥2M,, is multiplied in the first of Eqs. 18. Consequently the 7, curve is 
proportional to the torque, and the 6, curve to the unit angle of twist. 


be seen from Fig. 9 that the two curves greatly resemble each 


other. 
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case of a vanishing resistance to warping their ordinates would be proportional. 
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Fic. 10. 
Moment distribution 7p= —3Mtp//2 Pa, and unit angl 
twist distribution 6, (Eq. 39¢). 
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The calculation of the ‘* force ’’ terms (those containing Pa) is quite similar. 


From the end condition one obtains 
B= —.111 2/3) Pa? /2ET 


and thus one may write 


(B:!)p = — 9p 2Pa? /6EL 


(39a 
(39)) 


(39¢) 
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where 6, is a numerical factor. Its values are plotted against (s/a) in Fig. 10. 
The figure also contains calculated values of 7, which is defined as the factor 
by which — ¥ 2Pa/3 is multiplied in the first of Eqs. 18. In Fig. 10 the shape 
of the 6, and 7, curves differs considerably ; in the case of a vanishing resistance 
to warping their ordinates would be proportional. 
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Fic. 11. 
Plot of fy of Eq. 41 against g=GC/EI,, and h=GC/El). 


Principal axes of inertia as in Fig. 6. Effect of non-uniform 
warping considered. 


With the aid of the values of £;,' so obtained one can calculate M,,._ First the 
values of (@;')y and (8,')p must be multiplied by the corresponding values of t, 
in order to obtain the z-components of the rotations due to torsion. With these 
products plotted against s, a diagram is obtained in which the areas below the 
two curves are proportional to the angles of twist due to the ‘* moment ’’-terms, 
and the ‘* force ’’ terms, respectively. Carrying out the numerical calculations, 
one obtains 

a 
| =.0226y /2ED~.01275 2Put ‘CET (40) 
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remain the same as in Sec. 4. Since a?=1/K?=EI‘/GC one may write Eq. 40 
in the form 


a 
=.016 M,.a/GC—.003 Pa?/GC ; . (40a) 


The numerical factors on the right-hand side of Eq. goa are 16.6 times, and 
18.5 times smaller, respectively, than the corresponding terms in Eq. 19a. This 
shows that the torsional stiffness of the bar was increased to such an extent by 
the resistance to non-uniform warping. The expression for the redundant 
moment is 


where 
fu=(-045 + 1.259 — 4h) 
and g and h were defined in Section 4. Computed values of fy are shown in 
Fig. 11. The figure differs materially from Figs. 6, 7, and 8. 


g. THE EFFEcT oF NoN-UNIFORM WARPING UPON THE REDUNDANT MOMENT. 

If the space-curved bar is of thin-walled open section, its rigidity in simple 
torsion is very much smaller than that in bending. Thus in monocoque aircraft 
structures the usual values of g and i are around, or below 1/100. A comparison 
of Figs. 6 and 11 shows, however, little difference in the value of fy in this 
region. For very small values of g and h Eqs. 20a and 41a yield .208, and .188, 
respectively, for fy. It seems, therefore, that the value of the redundant moment 
is hardly affected by the increase in the rigidity of the bar due to the effects of 
non-uniform warping, although the calculations of the preceding article showed 
an increase of the rigidity to over 16 times its original value. 

The reason for this interesting phenomenon is that in the example the non- 
uniform warping had just about the same stiffening effect in the case of the 
external moment, as in that of M,,.. This is true of many practical problems, 
but not necessarily so. If the moment distribution due to the external loads 
differs materially from that due to the redundant moment, the stiffening effect 
may be different for different values of the ratio K7=GC/E; the value of the 
redundant moment will then also turn out differently. 


In order to get an insight into the numerical effects of different loadings upon 
the redundant moment, the bar shown in Fig. 12a is now investigated. The bar 
is rigidly fixed against torsion and warping at its left end a=L. At its right 
end, x=0, it is held fixed against rotation but it is allowed to warp freely. It is 
loaded with a distributed torque m expressed in in. Ibs. per in. according to the 
formula 

m=(u/sinhn) sinh (nr/L). (42) 
As the formula shows, the maximum value of m is u, and this value is reached 
at the left end. The shape of the distribution curve varies with the parameter 1 
and is shown in Fig. 12h for three different values of the parameter. 

If one chooses for the unknown the moment M, exerted upon the bar by the 
fixation at the right-hand side support, the ‘* simple *’ system is a cantilever bar 
with its left end rigidly fixed against rotation and warping. In this the bending 
moment due to the distributed torque is 


M=(pnL/n sinh n) [cosh (na/L)—1] ‘ (43) 
Substituting this value into Eq. 31a, one obtains 
sinh n) [cosh (nz/L)—1] (44) 


The end conditions are :— 


B' =o when 
B"=o when 
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The solution of Eq. 44 is 
=(uL/GCn sinh n) { [1 —(K?L? cosh n)/(K*L? — n*)] (cosh Kx) /(cosh KL) 
+ K*L? cosh (nx/L)/(K*L?—n?)—1 } (45) 
provided that lige) 
The relative angle of rotation 8, of the two end sections due to the distributed 
torque can be obtained by integrating the right-hand side of Eq. 45 between 
the limits =o and r=L. One obtains 
=(uL?/GCn sinh n) { [1 —(K?L? cosh n)/(K?L? — tanh KL 
+ K?L? (1 /n) sinh n/(K?L?—n*)—1 } (46) 
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Fic. 12. 
(a) Bar fixed at both ends against rotation, at left end against warping. 
(b) Distribution of applied torque of Eq. 42. 


In a similar way one may calculate the relative angle of rotation By of the 
two end sections of the bar due to a constant M,. The result is 
The redundant moment M, may now be calculated from the condition that the 
relative rotation of the two fixed ends must be zero :— 
Bu + By =o (48) 
Solving for M, one obtains 
M,=(nL/n sinh n) { [K?L? (1/n) sinh n | /(A?L*—n?) 
[1 —(K?L? cosh n)/(K?L?—n?)](1/KL) tanh KL } / { 1—(1/ KL) tanh KL } (49) 
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wl 


Next the value M,, of the redundant moment is calculated for the case when 
the bar does not offer any resistance to warping. The unit angle twist is then 
proportional to the moment :— 


B'=M/GC=(nL/GCn sinh n) [cosh (nz/L)—1] (50) 
The relative angle of twist 8, due to the external loads is 
B,.=(uL?/GCn sinh n) [(1/n) . (52 


The relative angle of twist 8,, due to the constant moment V,, is 


By (52) 
Applying Eq. 48 and solving for M,,, one obtains 
Moo= (ul sinh n) [1 —(1/n) sinh n | (53 


The effect of warping may now be calculated by determining the ratio M,/M,,. 
This was done for the three loading conditions represented by the values 1, 
and 10 of the parameter n. Computed values may be found in Fig. 13. 

Mo 
Moo 


8 


O é 4 6 8 10 
KL 


Fig. 13. 
Ratio of the unknown moment M, to the unknown moment 
M,., the latter calculated without consideration of the 
effect of non-uniform warping. 


Eq. 49 is not suited to numerical calculations when KL is very small. In such 
cases an expansion in power series yields the simpler formula 
M,=(pL/n sinh n) { 1 —(3/n*) [cosh n—(1/n) sinh n] } (54) 
It was mentioned above that Eq. 45 was not the solution of Eq. 44 when 
KL=n. In this case the solution is 
8'=(uL/GC sinh n) { (4) cosh Ka — (4) Ka sinh Kx —1 
— [{(4) cosh KL — (4) KL sinh KL —1] (cosh Kx)/(cosh KL) } (55) 
Integrating the right-hand side of Eq. 55 between the limits r=o and r=. 
and solving for M,, one obtains 
M,=(uL/n sinh n) { KL+2KL cosh KL —2 sinh KL 
— (4) sinh 2KL } / { 2KL cosh KL [1-(1/KL) tanh KL] } (5) 
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Fig. 13 contains computed values of M,/M,, in the region 
oSkKLS10 
The usual range in monocoque aircraft structures is within this region, as is 
shown in Sec. 11. The figure shows that there are considerable changes in the 
values of M,/M,, when the loading and AL are varied. 


10. CALCULATION OF THE UNIT STREss. 

After the statically indeterminate quantities have been determined with the 
aid of considerations of continuity, the calculation of the stresses in the cross 
sections does not present major difficulties. Stresses are caused by the direct 
force and the shearing forces as well as by the bending and twisting moments 
in each section, and the determination of the maximum combined stress is rather 
laborious. In this section the calculation of the stresses caused by torsion in 
thin-walled open sections is discussed in some detail. 

As shown in the Appendix (see Eq. A14, but substitute 8 for @), the direct 
stress o due to non-uniform warping is given by the equation 

o=EwB," 

If for 8,' an analytical expression is known, 8," is easily calculated by differentia- 
tion. If, however, 8,/ was obtained from a numerical, or graphical integration, 
the differentiation cannot be carried out analytically. In such a case two 
procedures are possible. The easier one consists of a graphical differentiation 
of the £,/=8,/ (s) function. The starting point for the more accurate one is 
Eq. 34. If one substitutes the expressions of Eq. 33a and Eq. 35 into Eq. 34, 
and then differentiates with respect to s, one obtains 


8," =K (A sinh Ks + B cosh Ks) — (1/2ET) fesfensa, (8) ds 


Since the values of the two integrals in Eq. 57 are known for a number of points 
along the centre line of the bar (that is for a number of values of the parameter s) 
from the numerical integration, and since the values of the constants A and b 
were also determined when £,' was calculated, the computation of 8,” from Eq. 57 
does not involve much additional work. The vaiue of the warping w can be 
calculated from Eq. A6. (See also Sec. A7.) 

In order to obtain numerical values for an appraisal of the importance of the 
effects of non-uniform warping upon the strength of the bar, the example of 
Secs. 4 and 8 is continued here. For the calculation of the stresses it is necessary 
to choose a definite section. The choice made here is a channel 2 ins. by 1 in. 
by .o4 in. Its sectional properties are as follows :— 

I, =.1067 in.‘ 
C= .0853 x 107° in.* 
in.-? 
I, =.0166 in.4 
I’ =.01166 
Assuming the length of one-quarter of the bar 
a =20 ins. 
one obtains 
Ka=1.06 
This value is close enough to unity so that the results of the calculations of 
Sec. 8 may be used without appreciable error. With 
g=GC/EI,=.308 x.10-* and 
h=GC/EI,=1.97 x 107° 
one obtains from Eq. 42a 
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and thus Eq. 41 becomes 
M,,=.175 Pa P (58) 


From Eq. 57 one can now calculate 8,” making use of the numerical values 
computed in Sec. 8. The results may be presented in the form 


-.0025 


7 -0020 -08 
-.0010 -.04 


.02 
-.0/0 


14. 
Numerical example of Sec. 10. C,=M,/Pa; 
C'=2GCB,'/ / 2Pa; C"=2GCB," 2P. 


In Fig. 14 values of C” are plotted against the parameter (s/a). The same 
figure contains two more curves. One of them represents the distribution of 
the torque M,=C,Pa. The values of its ordinates are the sums of the values 
of the ordinates of curve M,,/Pa, and of those of curve M,y/M,. of Fig. 5 
multiplied by .175. The other curve represents the distribution of the unit angle 
of twist B,/. It was calculated from the values of 6, and @, shown in Figs. 9 
and 10, respectively, with due consideration to Eq. 58. The ordinates are the 
values of C’ of the equation 

B/=C'J/2Pa/2GC . ‘ ; (60) 
With the aid of formule developed in Sec. AZ one gets the following values for 
the channel section used :— 


in. 

w =.625—8 sq. in. for the upper flange 

w= .1875 (28—4) sq. in. for the back (61) 
W = 3.375 —8 sq. in. for the lower flange 

d=.25 in. 
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The maximum value of w is reached at the free edges. It is 
Wmax = +-625 sq. in. 


Thus the maximum direct stress due to non-uniform warping is given by the — 


equation 
= Ib. per sq. in. 
where x is a function of the parameter s. Values of o,/P are plotted against , 
(s/a) in Fig. 15. ] 
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FIG. 15. accu 


Numerical example of Sec. 10. Direct stress in edge fibre. 
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For the sake of comparison the direct stress caused by bending is also computed. 
The resultant bending moments can be obtained from the curves of Fig. 5 by 
adding up the values of the moments with the subscript P to the corresponding 
ones with the subscript M, the latter multiplied by .175 on account of Eq. 58. 
Dividing the moments so obtained by the section moduli 

Z,=.1067 in.* and Z,=.0221 
one obtains the direct stresses o, and o,, due to the bending about the two 
principal axes of inertia. These are also shown in Fig. 15. 

Each of the three curves o,/P, o,/P, and o,/P represents the variation of the 
direct stress along one free edge of the channel caused by torsion with non- 
uniform warping, by bending about the principal normal, and by bending about 
the binormal of the median line of the bar, respectively, assuming the directions 
of the principal normal and of the binormal to be the directions of the two 
principal axes of inertia. It may be seen that the direct stresses caused by 
torsion and by bending about the principal normal are approximately of the same 
magnitude while those caused by bending about the binormal are about ten times 
higher. If one would turn the channel through a right angle so that the greater 
moment of inertia would correspond to the binormal, and the smailer to the 
principal normal, the stresses would naturally be reduced. Such a change in the 
orientation of the channel would couse only a slight change in the value of fy, 
and Eq. 58 would become 

M,,=-185 Pa 

Thus the values of the bending moments, and of 8,”, would not change materially. 
Since, however, the ratio of the section moduli is 4.82, the stresses o, caused 
by bending about the principal normal would increase roughly to five fold, those 
caused by bending about the binormal would decrease roughly to one-fifth their 
original values. As o;, would remain unchanged, its maximum value would be 
roughly one-fifth of the maximum bending stress. This is certainly high enough 
not to disregard except in rough approximate calculations. 

The shearing stress may be calculated from Eq. A24 of the Appendix. Writing 
8 tor @, one has 


s 


7= ds 
The integral can be calculated if one substitutes the value of w from Eqs. 61 :— 
ds=.625 s* for upper flange 


ds =.1875 8° —.75 s+.0875 for back (62) 


G ds = 3.375 8—.5 s°—5.5 for lower flange 


The value of this integral is proportional to the shearing stress, the distribution 
of which is shown in Fig. 16. The figure also contains the direct stress 
distribution caused by non-uniform warping. 
The maximum value of the integral is .1953 in.* at a distance of .625 in. from 
the free edge. Hence the value of the maximum shearing stress is 
It is possible to calculate £,/’ by differentiating, analytically or graphically, 
the function 8,” ='8," (s) previously obtained and shown in Fig. 14. It is more 
accurate, however, to calculate it from Eq. 31a. Solving this for 8,/’, one obtains 


= K2B/ (5 /ET) M, (8) 


5 
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Since the values of both 8,’ and M,(s) were previously computed for a number 
of values of the parameter s and plotted in Fig. 14, ,!” is easily obtained for 
the same values of s. The maximum shearing stress can then be calculated from 
Eq. 63. 

At the same time shearing stresses are caused by the part of the torque that 
is resisted by simple (Saint-Venant) torsion. It is known from the theory of 
elasticity that the shearing stresses in a thin-walled channel section are very 
nearly the same as those in a narrow rectangular section the height of which is 
equal to the developed length of the median line of the channel, and the width 
of which is equal to the thickness of the channel. The maximum shearing stress 
occurs in the middle of the longer side of the rectangle. The value of the shearing 
stress due to simple torsion at s=.625 in. and n=t/2 does not differ much from 
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Fic. 16. 
Stress distribution in a channel section due to non-uniform warping. 


this maximum, as may be concluded from the statements made in Sec. A2 and 
from Fig. A3. Assuming that it is the same, and calculating the maximum stress 
for a rectangular section of a height h=2 ins. and a width t=.04 in., one may 
write 


7=£8,'/tG@ =6630 C'P Ibs. per sq. in. : (64) 
where C’ is the non-dimensional coefficient plotted in Fig. 14. 

The variation of the shearing stress in the most stressed fibre (in shear) is 
shown in Fig. 17. The total shearing stress is the sum of the two components. 
It is seen that in the case under consideration the maximum shear associated 
with non-uniform warping is roughly equal to that associated with simple (Saint- 
Venant) torsion. The maximum absolute value of the sum of the two is about 
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37 P lbs. per sq. in. If one calculates the maximum reduced stress for failure 
according to the Mohr-Guest theory, 


Trea = (o? 
one finds that the shearing stress has a considerable influence upon failure. 


Trax 
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17. 
Numerical example of Sec. 10. Maximum shearing stress. 


For a check, the torque caused by the two different kinds of shearing stress 
is calculated for (s/a)=.8. The shearing stress in simple torsion is —15 P lbs. 
per sq. in. The corresponding torque is 

My = (1/3) = 2.139 — 032 P in, Ih. 
The maximum shearing stress associated with non-uniform warping is — 20.4 P Ibs. 
per sq. in. Since the average shearing stress in the flange is .75 times the 
maximum, as may be calculated from the first of Eqs. 62, the total shear force 
S, in the flange is 
S,= —.75§ X 20.4 1 x .04=.612 P 
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An equal and opposite shear force is acting in the lower flange; the moment of 
the two 

M,= —1.224 P in. lb. 
The total moment is the sum of the two 

tot = — 1.256 P in. Ib. 
Here the moment due to simple torsion is about 2.5 per cent. of the total; 
nowhere does it exceed 7 per cent. in this example. Nevertheless, as is seen 
from Fig. 17, the stresses associated with non-uniform warping are about equal 
to those due to simple torsion. 

The torque at (s/a)=.8 is, from Fig. 14, 

M,= —.0623 Pa= —1.246 P 
The deviation from the formerly calculated value is less than 1 per cent. This 
is regarded as satisfactory, since the calculations were carried out with the 
slide-rule. 

It should be mentioned that in the example the numerical values obtained for 
the shearing stress associated with non-uniform warping are so high that the 
assumptions underlying the theory presented in the Appendix cannot be considered 
as satisfied. Thus the values of the stresses calculated for the channel section 
chosen are not reliable. Fortunately the sections actually used in aircraft con- 
struction are generally heavier than the section in the example; in them the 
shearing stresses associated with non-uniform warping are smaller, and conse- 
quently the theory of the Appendix applies. 


11. SIMPLIFIED PROCEDURE. 


A. If the shape of the centre line of the curved bar is reasonably similar to the 
intersection of two straight circular cylinders, then Marguerre’s diagrams 
may be used for the determination of the unknown moment M,,, of the 
moment distribution, and of the maximum moment, provided that the 
orientation of the principal axes of inertia of the section is such as assumed 
by Marguerre, or if not that it has little effect upon the moment distribution, 
and provided that the effect of non-uniform warping upon the unknown 
moment M,, is negligibly small. 

a. The directions of the principal axes of inertia permit the use of 

Marguerre’s diagrams :— 

1. If one of the principal axes of inertia is always horizontal, as assumed 
by Marguerre. This is the case, of course, if the section is circular, 
or circular cylindrical, since then every axis is a principal axis; 

2. Or if the ratios g=GC/EI,, and h=GC/EI,, (where the numerator 
is the Saint-Venant torsional rigidity, and the denominators are the 
bending rigidities about the two principal axes of inertia) are both 
very small, say around or below .or. 

For the solid circular, and the circular ring section these ratios 
are .77. For a solid, narrow rectangular section of a width b and 
a height h one has 

C=(1/3) b*h I,,=(1/12) bh? 1,,=(1/12) 

g=-385x4(b/h)? h=.385 x 4=1.54 
In the case of common thin-walled open sections g and ht are usually 
small. For instance, in the case of the angle section shown in 
Fig. A5 the moment of inertia about the horizontal principal axis is 
T=(1/12) [1 +6 (b/h)] 
The torsional constant C is given by the formula 


C=(1/3) ht® [1+2 (b/h)] 
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Thus 
g=1.54 (t/h)* [1 +2 (b/h)]/[1+6 (b/h)] 
For the heavy channel 1 in. by } in. by } in. g becomes .o12; for 
the light channel 2 ins. by 1 in. by yg in. the vaiue of g is .00075. 
Although the value of h is greater, its order of magnitude is the 
same, 
b. Marguerre’s diagrams should not be used because of the effects of non- 
uniform warping if the resistance of the section to non-uniform warping 
is considerably greater than its resistance to Saint-Venant torsion, and :— 
1. If there is reason to believe that non-uniform warping influences the 
distortions caused by the external moments differently from the way 
it influences the distortions caused by the unknown moment (see 
Sec 9); 

2. Or if g and ht are not very small. 

It may be seen from Secs. 7 and 8 that the effect of non-uniform 
warping upon the rigidity of a bar depends upon the value of the 
non-dimensional coefficient Ka. The probable range of this co- 
efficient is, therefore, now established. In the case of the channel 
section of Fig. As the warping constant I’ is given in Eqs. A47 and 
A47a. Hence 

(GC /ET)=1.54 (#?/h4) (h/b)® [1 +2 (b/h)] [1+ 6 (b/h)]/[2+ 3 (b/h)] 
The value of K=(GC/EI’)} is, therefore, .66 for the above-mentioned 
‘‘ heavy ’’ section, and .083 for the ‘‘ light ’’ section. Taking the 
length a of one-quarter of the centre line of the bar 15 ins. in the 
former case, and 12 ins. in the latter, one obtains the approximate 
limits of Ka :— 
1X KaS10 
In the case of the angle section shown in Fig. A6b one may 
calculate the value of K with the aid of Eq. Az7o:— 
K=2.15 (a+ + 

Thus for most angle sections used in aircraft K is greater than 
unity, and for most space curved rings Ka is greater than ten. 

Numerical calculations show that the increase in stiffness of a 
twisted bar due to its resistance to non-uniform warping is 
negligibly small under normal loading conditions if the value of Ka 
is above 10; as was found in Sec. 8, it is very great if Ku is around 
1. In the range from one to ten the effect of non-uniform warping 
upon the stiffness changes very rapidly. 


. The procedure developed in Secs. 3 and/or 4 should be used for determining 


M,, and the moment distribution, if Marguerre’s diagram cannot be used 
because of the shape of the centre line of the bar, or because of the reasons 


enumerated under A, and if it can be assumed that the effect of non-uniform 
warping on M,, is small. 


The procedure developed in Secs. 7 and 8 should be used if the effect of non- 


uniform warping upon M,, is great. ‘This is the case when the conditions 
stated under B hold. 


. For an accurate determination of both the direct and shearing stresses the 


distribution of the angle of twist must be calculated from the differential 
equation given in Sec. 7. This done, the stresses may be calculated as 
shown in Sec. 10. For an approximate calculation, the effects of non- 
uniform warping may be neglected if Ka is great; if it is very small, and 
thus the shearing stresses associated with non-uniform warping are very 
large, the theory developed in the Appendix cannot be relied on. 


¢. In the calculation of the distortions, as a rule, the effects of non-uniform 


warping must not be neglected if a thin-walled open section is used. 
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APPENDIX. 


A1. SOLUTIONS OF THE PROBLEM OF TORSION. 


Since de Saint-Venant established the basis of the present day theory of the 
torsion of straight bars, the shearing stress distribution in and the angle of twist 
of twisted bars has been calculated for a great number of different cross sections. 
Some important assumptions underlying these solutions are that the twisting 
moment is constant along the bar, it is applied by means of two equal and 
opposite couples at the two end sections of the bar, that these external couples 
are the resultants of shear stresses distributed over the end sections in the same 
way as are the shear stresses in any and every cross section of the twisted bar, 
and that the end sections are free to warp. This last assumption means that the 
points in the plane of the end sections are not constrained to remain in the plane 
when the bar is twisted. In this Appendix stresses, angle of twist, etc., cal- 
culated in such a way will be marked ‘‘ in torsion according to Saint-Venant,’’ 
or ‘* according to Saint-Venant.”’ 


Fig. Ar. 
Torsion with non-uniform warping. 


If these assumptions do not apply, the stress distribution in and the torsional 
rigidity of the twisted bar may be greatly affected. A problem of this kind was 
first discussed in 1g05 in Russian, and in 1910 in German in papers written by 
S. Timoshenko,* in which he calculated the effect upon a twisted I-beam of 
preventing its fixed end section from warping. Fig. A1 shows that in such a 
case the distortions of the twisted bar involve a so-called differential bending, 
namely, the bending of the two flanges in opposite directions, in addition to the 
twisting according to Saint-Venant. Consequently direct stresses as well as 
shearing stresses are caused by the torque. This more complex problem is 
governed by the differential equation 

T =GCq! — (h?/2) ‘ (At) 
where T is the torque, GC the torsional rigidity of the bar according to Saint- 
Venant, E Young’s modulus, J, the moment of inertia of the cross section of 
one flange about an axis through its centre of gravity parallel to the y-direction, 
@ the angle of twist, and the dash denotes a differentiation with respect to a. 


Eq. Ar was developed in Timoshenko’s above-mentioned papers. A derivation 
and discussion of this equation will not be given here, since they can be found 


in Ref. 5. 


* See Ref. 4. Also p. 282 of Ref. 5. 
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A2. STRESSES IN A Bar oF THIN-WALLED OPEN SECTION TWISTED ABOUT ITS 
CENTROIDAL AXIS. 


It is not difficult to establish an equation similar to Eq. A1 for a twisted bar 
which has a thin-walled open section of the kind used in aircraft construction. 
Two papers written in German by H. Wagner and R. Kappus, respectively, and 
published by the N.A.C.A. in English translation, deal with this subject. (See 
Refs. 7 and 2.) Notwithstanding the fact that the English translations are 
easily available, a derivation of the differential equation of torsion with non- 
uniform warping of thin-walled open sections will be given here since the publica- 
tions are not easy to follow. One cause of this lack of clarity is that reference 
is made to a paper by W. Liicker for some of the proofs, and this paper is 
unknown to the writer and not available in this country. 


TANGENT ATS 


Fic. A2. 


Median line of a thin-walled open section. 


Fig. A2 represents the median line of the cross section of a thin-walled open 
section bar. The z-axis is coincident with the longitudinal axis of the bar and 
points toward the reader so that the axes 2, y, 2 represent a right-hand system. 
The position of a point of the cross section is essentially determined by its 
distance s from the point s=o, measured along the median line, since the wall 
thickness t is considered small. The positive direction of the tangent is deter- 
mined by the rule that the tangent should point in the direction of increasing 
values of the co-ordinate s. The perpendicular distance r, of the centroid from 
the tangent is taken positive if the moment of the tangent vector about the 
centroid is as shown in the figure, that is when it corresponds to a rotation from 
the y-axis toward the z-axis. Furthermore, the bar, the cross section of which 
is shown in Fig, A2, is assumed to be straight, of constant cross section, and 
made of homogeneous isotropic material which follows Hooke’s law. 


Let it be assumed that two cross sections of the bar a distance dz apart rotate 
about the a-axis, which passes through the centroid, by an angle d@ relative to 
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each other. Such a rotation will cause shearing stresses to appear in the cross 
sections which may preferably be calculated with the aid of the membrane analogy.* 
Fig. A3 shows the contour lines of «a membrane, attached to the edges of an 
open section which has a wall thickness t, as they develop under uniform lateral 
pressure. These contour lines can also be considered the lines of shearing stress 
of the twisted bar. It may be seen that no shearing stresses are acting along 
the median line of the thin wall. From this, however, it follows that the 
originally right angle between any line element ds of the median line and the 
direction of the longitudinal fibres of the bar, originally going in the z-direction, 
must remain a right angle after distortions, since by the laws of the theory of 
elasticity the distortion of this angle is proportional to the shearing stress which 
acts in the cross section in the direction of the median line. This requirement 
permits the calculation of the warping of the cross section. 


MeEDIAN LINE 


Fig. A3. 


Lines of shearing stress in a twisted thin-walled open section. 


Let u denote the displacement of a point of the median line of the cross section 
from its original position in the x-direction, and v that in a direction tangential 
to the median line. Then the change y of the originally right angle between 
the longitudinal fibres of the bar and the tangent to the median line of the cross 
section is given by the known equation of the theory of elasticity 

y= (dv /dx) + (du /ds) (A.2) 
This change, however, must vanish as stated above and therefore u can_ be 
obtained by an integration 


us — | (00/02) ds-+F (a) 


if dv /dr is known. In Eq. A.3 F (x) is an arbitrary function of x. Since v is 
the displacement of a point of the median line in the direction of the tangent 
to the median line, and since the relative displacement of the two neighbouring 
sections, by assumption, consists of a rotation about the centroid of the section, 


* See p. 239 of Ref. 6. 


wl 


= — | 
/ 
\\ 
\\\ \ \ \ 
\\ \ \ \ \ \ 
| 
| | | j 
| 
/ / / 
yy, | 
/ /, / 
7 { 
| He 
of 
int 


STRESSES IN SPACE-CURVED RINGS IN MONOCOQUE FUSELAGES. 67 


where the dash again denotes a differentiation with respect to 2, and ry is a 
function of the co-ordinate s alone. When dv is introduced from Eq. A.4 into 
Eq. A.3 the longitudinal displacement wu becomes 


u= of res +F (x) (A.5) 
where F (x) is a function of x alone. 

In further calculations it will be convenient to measure the longitudinal 
displacements from an average-displacement-plane rather than from the original 
position of the points, and to define a quantity wy, as such a displacement 
corresponding to a unit rate of change of the angle of twist :— 


Wo= [reds +C . . (A.5) 


Sh 


Fic. A4. 
Graphical determination of the warping (Eqs. A.0). 
where 
L 
[red 


Here the warping we is a function of s alone and is determined by the shape 
of the median line of the cross section. It can be readily obtained by a graphical 
integration, since {r,ds is the double area enclosed by the segment of the median 


* Throughout this paper this form of setting the double integral is the same as r.dsds 
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line between o and s, and by the two straight lines connecting these two end 
points to the centroid of the section. (See Fig. A4.) 
The strain in the «x-direction is, by definition, 
= (du /dz) 

and the longitudinal stress 

o,=E (du/dz) (A.7) 
Eqs. A.5, A.6, and A.7 yield 

(x) . : (A.8) 
where f (x) is an arbitrary function of 2 replacing the more complex expression 
E [F' (x)—Co"]. It can be determined from the requirement that for the case 
of pure torsion the resultant axial force must vanish. If one considers that the 
stress must be substantially constant across the small constant wall thickness ft, 
the resultant axial force X becomes 

L L 


X= | o,ds [weds +tLf (zx) 


Since from Eq. A.6 it follows that 
L 


the condition of a vanishing axial force X is 
f (x)=0. 
Hence the axial stress distribution is defined by the equation 


where w, is the quantity defined in Eq. A.6. 


A3. THe CENTRE oF TWIST. 

The stress distribution calculated above corresponds to a vanishing resultant 
axial force X, but in general the resultant moment of such stresses does not 
vanish about either the y- or the z-axis. The equilibrium conditions cannot be 
satisfied, therefore, if the external force system consists of a pure torque. This 
contradiction may be eliminated by assuming that the bar does not twist about 
an axis through the centroid but rather about an axis through a point T. The 
co-ordinates yy, 2, of this point T can be determined from the requirements that 
the resultant moment should vanish about both the y- and z-axis. 

For a rotation about an axis through T the distance r, from the centroid must 
be replaced by the distance r, from point T in the formulas developed in Sec. A2 
(see Fig. Az) :— 

and a warping w, may be defined in analogy to w.:— 
L 


s 
| + (1/L) (as [rts 
: 


8 
=—[ (re- 4—Z_ Sin x) ds (A.12) 
0 


8 
+(1 (ro — COS Z_ Sin x) ds 


L 


|_| 
ts 
{ 
ig 
5 a 
3 
tl 
4 
| | wl 
mc 


A .12) 


STRESSES IN SPACE-CURVED RINGS IN MONOCOQUE FUSELAGES. 69 


Considerations similar to those made under Sec. A2 lead to the result that the 
axial stress o, is 
and that with such a stress the force vanishes. 
With the aid of the identities 
ds cos 2a=dz and dssinx=—dy : (A.15) 
Eq. A.12 can be transformed :— 


s 


— [res + Yq | dz — zy | dy 


0 


Ls z 
+(1 [| ret de [ay] ds 
e 


L 
= We t+ Yq (2 — 20) — (Y— Yo) — (1 /L [Yap — — (Y— Yo) ds 


since the sum of the original first and fourth terms is w, on account of Eqs. A.6. 
This equation can be further simplified by evaluating the definite integrals :— 
L 


(1/ 

0 
This is true since z and y are measured from centroidal axes and since the 
symbols with a subscript are constants. With these values the warping wy, 
caused by a rotation about a longitudinal axis through point T becomes :— 

Wy =Wo + — ‘ (A.16) 
which equation means that the warping wy, differs from the warping w, by an 
amount corresponding to a rotation of the plane of the cross section an axis 
contained in the plane. Upon substitution into Eq. A.14 of wy from Eq. A.16 
the stress becomes :— 

The unknown co-ordinates y, and z, can be calculated Sai the requirement 
that the resultant moment must vanish about both the y- and z-axis :— 


L L 
Oo 


If from Eq. A.14 ©, is introduced into these ie ns they become :— 


where w, is the warping given by Eq. A.16. With the symbols for the second 


moments of the cross section 
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warping moments 


and for the 


and Fat | ‘ (A. 20) 
Eqs. A.18 become :— 
—F, 


Eqs. A.21 can be solved for y, and :— 
(Falyy — yy — Lye”) (A.22) 

The point defined by the co-ordinates yy, 2, given by Eqs. A.22 is called the 
centre of twist. The cross section of the twisted bar must rotate about an axis 
passing through the centre of twist if the external load is a pure torque and it 
the bar cannot warp uniformly so that in addition to the shearing stresses axial 
stresses are caused in the cross section. It should be noted that in Eqs. A.19- 
A.22 y and z are measured from the centroidal axes, and uw, must be calculated 
from Eqs. A.6 for a rotation about the centre of gravity of the section. 

In the above developments it was tacitly assumed that the shearing stresses 
caused by non-uniform warping were small as compared to those due to Saint- 
Venant torsion. Only in this case is Fig. A3 representative of the shearing stress 
distribution in the twisted bar. 


A4. SHEAR STRESSES DvuE TO RESTRAINED WARPING. 
An equilibrium condition of an infinitesimal element of the thin wall is 
(0c, /0x)+(07/ds)=o (A.23) 
This differential equation permits the calculation ot the shearing stresses in the 
thin wall, since the distribution of the axial stresses is already known. By 
introducing o, from Eq. A.14 into Eq. A.23, and solving the equation, the 
following expression for the shearing stress is obtained :— 
s 
T= fire, (Ox) ds +y (x) 


= +g (x) 

0 
where g(x) is an arbitrary tunction of «2, to be determined from the boundary 
conditions. Since no external shearing forces are introduced along the edges 
s=o and s=L, the shearing stress must be zero at the edges :— 

7=o for s=o and s=L. 
These conditions can be satisfied if 
g (x)=0, 

since Eq. A.13 holds true. Consequently the shear distribution is given by 


s 


0 


It can easily be shown that the resultant of these shearing stresses is a pure 
torque. The y-component of the resultant is 
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This integral can be evaluated by integrating by parts :— 


L 
- 
{ [yf |+ | wryds }. 
0 0 


In this equation the upper limit of the first term is zero in accordance with 
Eq. A.13, and the lower limit is equal to zero for obvious reasons. The second 
term vanishes in accordance with the first of Eqs. \.180. Hence 


: (A.25) 
Similarly 
L 
cos A cos ads = — cos | wyds 
A 
becomes 


L 


L 
S,=—Eo'"t { wads] } 
oO 


which vanishes for reasons similar to those shown above, giving 


Having thus ascertained that there is no resultant shear force due to restrained 
warping, one can proceed to the calculation of the resultant torque : 


L L s 
— Eq!"t | wis 
o 
Integrating by parts, one obtains the integral 
I. 
L 
T=—Eol"t { [- + } 
) 


8 
Since {wyds vanishes for the upper limit by reason of Eq. A.13, and for the 


lower limit for obvious reasons, the above equation can be rewritten in the form 
where L 


0 
I’ will be called the warping constant.* From the definition of w, it follows 
that the warping constant is a cross-sectional constant, it has the dimension 
of a length to the sixth power, and is always positive. 


As. THe DIFFERENTIAL Equation or Torsion with NoN-UNIFORM WARPING FOR 
THIN-WALLED OPEN SECTIONS. 


The findings of the above sections show that in general if the torque acting 
upon a straight thin-walled open section bar varies along the bar, or if the free 
warping of any cross section is prevented, axial stresses will appear in the thin 
wall in addition to the shearing stresses according to Saint-Venant. These axial 
stresses will always be accompanied by shearing stresses, the distribution and 
magnitude of which are entirely different from those according to Saint-Venant. 
The shearing stresses associated with restrained warping add up to a pure 
torque which together with the resultant torque of the Saint-Venant shearing 
stresses balances the external torque applied to the bar. 


* This symbol was used by Goodier in Ref. 1. 
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Since the torque due to the shearing stresses according to Saint-Venant is 

T sy =GCq! (A.29) 

where GC is the torsional rigidity according to Saint-Venant, the connection of 

the external torque with the internal torque (due to stresses) can be written, 
with the aid of Eqs. A.27-A.20, 

This is the differential equation of torsion with non-uniform warping for thin- 

walled open sections and is very similar to Eq. A.1. If the geometry and the 

material properties of the bar, and the variation of the external moment are 

known, the angle of twist can be calculated by solving this differential equation. 


A6. THE SHEAR CENTRE. 


For numerical computations it is helpful to know that the centre of twist of a 
section is identical with its shear centre. This follows from the reciprocal 
theorem. It is interesting, however, to prove formally that the requirements 
for the shear centre of a thin-walled open section lead to the same formule as 
those found previously when the centre of twist was determined. The shear 
centre is defined in the theory of bending of non-symmetrical sections as the 
point through which an external shear force must pass if it is to cause bending 
alone and no torsion. Alternatively it can be stated that the resultant shear 
force always passes through the shear centre, if the section is only bent and 
not twisted. This second definition may be used advantageously to calculate the 
location of the shear centre. 

If the direct stress is known, the shear stress distribution can be calculated 
by solving the differential equation 

From this 
s 
T= | ds + G (x) 
where G (x) is an arbitrary function of z alone. 

The axial stress in a non-symmetrical beam bent by two moments M, and M, is 

given by the formula 

(1,,M, 4 190") 

—y + I5zMy) 1 ys”) 

where y and z are any two mutually perpendicular centroidal axes in the plane 
of the cross section, M, and M, the applied bending moments whose vectors are 
parallel to the y- and z-axes, respectively, and point in their positive directions 
(whereby the sense of the vectors corresponds to the right-hand screw rule), 
and the cross-sectional properties of the bar are defined by the integrals 


yeds 


Since 


(dM,/dx)=S, and (dM,/dz)=—S, . (A.31) 
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where S, is the resultant shear force acting parallel to the z-axis, and S, the 
resultant shear force acting parallel to the y-axis, the differential coefficient of 
the axial stress is 
(00, /Ox) =z / — Is) 
+y — / 12”) 
For a given bar and given loading every symbol but y and z of the above 
formula represents a constant. With the introduction of the symbols 


C,= (IyyS, — / — (A.32) 

the expression for the shearing stress becomes 
Cy) ds+G (a), ‘ (A.33) 


The boundary conditions are 


7=o for s=o and s=L. 
They are satisfied if 
G (x)=o0, 
since 


A 


1 
=o 
v A 
for the reason that y and z are measured from centroidal axes. Consequently 
Eq. A.33 reduces to 


Ss 


| Cy) ds ‘ (A.34) 


wherein C, and C, are constants defined by Eqs. A.32. 

Shearing stresses distributed over the cross section according to Eq. A.34 
add up, in general, to a resultant shearing force having a y- and a z-component, 
and to a resultant torque. With a view to checking Eq. A.34 the y-component 
of the resultant shearing force will now be calculated. 


L 
=— sin xdA= — sin ads 
A A 

L 
=| sin ads | +C,y) ds. 
0 
Integrating by parts, one obtains 


L 
L 
| { [y7]+ (C,2+C,y) ds}. 
A 
On the right-hand side the first term vanishes since 7 vanishes for both the 
upper and lower limit. The second term becomes 


C Ty, + Cylyy. 
With the substitution of C, and C, from Eqs. A.32 the integral becomes 
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In a similar way it can be shown that 


A 
Consequently the shearing stresses in the section, as given by Eq. A.34, actually 
add up to the total applied shearing force. Thus the only remaining task is to 
determine point S in the plane of the cross section about which the resultant 
torque T of all the shearing stresses vanishes. By definition this point is the 
shear centre of the cross section, since the moment of the shearing stresses 
vanishes identically only for the point through which the resultant shear force 
passes for any arbitrary choice of S, and S,. Denoting by r, the perpendicular 
distance of the shear centre from the tangent to a point of the median line of the 
thin wall of the section, one obtains the equation of the moment about the shear 
centre :— 


L L 


Integrating by parts, one obtains :— 
i, 


L 
T=-—t { [wz] + [ice +C,y) }. 


In this equation the first term on the right-hand side vanishes since 7 vanishes 
for both limits. Since C, and C, are entirely arbitrary for the reason that S, 
and S, may have any arbitrary value, the torque T vanishes identically only if 


L L 
| | Cayre.ds=o (A.37) 


These two conditions determine the location of the shear centre. Eqs. A.37, 
however, identical with Eqs. A.18a except for a constant factor. Consequently 
the requirements for the shear centre and those for the centre of twist are 
identical, and therefore the centre of twist must coincide with the shear centre. 


Az. EXAMPLE FOR THE CALCULATION OF THE WARPING CONSTANT. 


The first step in calculating the warping constant |’ is to determine the centre 
of twist. This may be done with the aid of Eqs. A.22 but in most cases it is 
more convenient to start from the definition of the shear centre. This definition 
was given and applied to the general thin-walled section in Sec. A6. Practical 
calculations for commonly used sections were carried out by Timoshenko in 
Art. 8 of Ref. 5. 

The second step is to determine the warping wy, for a rotation about the centre 
of twist. This can always be done easily for thin-walled sections of any 
arbitrary shape by the graphical method explained in Sec. A2 and shown in 
Fig. A4. When this is accomplished, a numerical or graphical integration of 
the integrand in Eq. A.28 will not present any difficulties. It should be noted, 
however, that the formule in Sec. Az, and Fig. Aq refer to a rotation about 
the centroid, and not about the centre of twist as required here. 

In the case of sections of a shape which can be represented by simple 
mathematical expressions, the warping constant I’ can be found analytically. 
This will be shown in the example of the channel section of Fig. As. On the 
assumption of a small wall thickness t the moment of inertia of the section is 


= (t/12) (6bh? +h?) i (A.38) 
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The distance e of the centre of twist from the back of the channel is 


according to the formula on p. 53 of Ref. 5. This can also be written 


Since ry remains constant along the two flanges and along the back, the first 
integral on the right-hand side of Eq. A.12 becomes :— 


r,ds= —(h/2)s when oSs Sb 


rpds=es —(b/2)(2e+h) when bSs< beh (A. 40) 


reds = —(h/2) [s—(2e+h)] when b4hSs5<2b+h 
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Fic. As. 
Warping constant of the channel section. 
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With these values the second integral on the right-hand side of Eq. A.12 
becomes 


L 8 
ds rls = (1/2) (A.41) 


The warping w, for a rotation about the centre of twist, therefore, is given 
by the expressions 


Wp=(h/2) (D-—s) when oSsSbh 
(e/ 2)(2s—L) when bSs=b)+h (A.42) 
Wp=(h/2)(L—D—s) when b+hSsS 
where L is the total length of the centre line of the section 
L=2b+h ‘ (A.43) 
and D denotes the expression 


The calculation of I’ according to Eq. A.28 is now straightforward. By 
symmetry 


L b bt+h 


The integration gives 


= /6) — 3Db? + b? | 


7 ds =(t/12) e*h*® 


After some algebraic iciicaeaibia the warping constant becomes 
=(th?/12) +2b*—6eb (b—e)] . (A.45) 
If the value of e is introduced into this formula from Eq. A.39a, the warping 
constant can be written as 
(th*b*/24) [1+ 3h/(6b+h)] (A.46) 
If this expression for I’ is introduced into Eq. A.27, it can be easily verified 
that the result is identical with Eq. (y) on p. 292 of Ref. 5. A more convenient 
form for numerical computations is 
where 
k=(1/12) [3 (b/h)*+2 (b/h)§]/[6(b/h)+1] 
The value of k is plotted in the chart of Fig. A.5. 


A8. THe ErFect oF A NoN-VANISHING WALL ‘THICKNESS. 


A refinement of the foregoing calculations may be accomplished by a considera- 
tion of the effect of a non-vanishing wall thickness. Wagner showed this in 
Ref. 7. This effect is of importance with sections for which I’ (Eq. A.28) 
vanishes, as is the case for instance with angle- and T-sections. 

The distance of a point of the section from the median line, measured 
perpendicularly to the tangent to the median line, is denoted n. The positive 
direction of the perpendicular may be obtained by rotating a vector which points 
in the positive direction of the tangent clockwise through a right angle; it is, 
therefore, the ‘* outward ’’ direction in Fig. Az. The distance of the perpen- 
dicular from the centroid is denoted r,,,; it is considered positive if a force pointing 
in the positive direction of the perpendicular causes a counter-clockwise rotation 
about the centroid. Consequently r,, is positive in the case of point s shown 
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in Fig. Az. The displacement of a point from its original position during the 
elastic distortions measured in the positive direction of the perpendicular is 
denoted by the letter p. 

From the basic equations of the theory of elasticity the change y of the 


originally right angle between the perpendicular to the median line and the 
longitudinal fibres of the twisted bar is 


= (dp /dx) + (Ow /dn) (A. 48) 
On the other hand it may me seen from Fig. A3 that a for a very limited 
region at the ends of the thin-walled section the direction of the Saint-Venant 
shearing stresses is parallel to the median line. If one assumes that torsion with 
non-uniform warping of a section with a non-vanishing wall thickness does not 
cause shearing stresses of any appreciable magnitude in the direction of the 
perpendicular, then the angle y=7/G must vanish. Hence 

n 


u= | dn+u, . (A.49) 
0 
where u, is the longitudinal displacement of a point of the median line. From 
the assumptions concerning the distortions of the bar, stated in Sec. A2, it 
follows that 


dp dx . (A.50) 
Substituting into Eq. A.49, one obtains 
n 
oO 


since /y¢ IS a constant when one integrates with respect to wu. 


The displacement u, of the median line was calculated in Sec. 2. It is given 
in Eq. A.5 in which it is represented by the symbol u without the subscript o. 
Substituting its value into Eq. A.51, and using an asterisk to denote that the 
displacement u (or any other quantity in future calculations) was determined 
with due consideration of the effect of a non-vanishing wall thickness, one obtains 

u*= ds — 9! Mot + F (x) ‘ (A.52) 

Since n is measured from the median line, it follows from Eq. A.51 that the 
average displacement across the wall thickness is the displacement at the median 
line. Consequently the warping w,* differs from the warping w, defined in 
Sec. 2 only by a quantity corresponding to the first term on the right-hand side 
of Eq. A.51 :— 

Wo" . : ‘ 
The warping w, is given in Eqs. A.6. The direct stress o,* can be calculated 
in the same way as in Sec. 2, and since the average stress across the wall is the 
stress at the median line, one obtains an equation analogous to Eq. A.10 :— 

If the cross sections of the bar rotate about the centre of twist 7' rather than 
about the centroid C, one may write 

We . : (A.5s5a) 
and 


o,*=¢"w,*E . ; (A.g5b 
The connection between the stresses corresponding to a rotation about the centre 
of twist and about the centroid can be established with the aid of Eq. A.16. 
Eq. A.55a may be written in the form 
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Since the y- and z-components of n are 

ny=ncos xz and n,=Nn sin x 
and further since 

Top — Yr SIN 44+ COS 
one obtains, with consideration of Eq. A.53, 


=U + (2 +z) — Zp (Y + Ny) 


cross section. 


Moment equilibrium requires that 


t/2 
fos (y+ny) dn=o 
—t/2 
and 
t/2 
|as (2+n,)dn= 
-t/2 


In this equation (y+n,) and (z+n,) are the co-ordinates of any 


point 


where y and z are co-ordinates of the median line measured from the centroid. 


(A.57 


of the 


(.A.58) 


where o,* has the value given in Eq. A.55b. These requirements permit the 


calculation of the co-ordinates y,* and z,* of the centre of twist :— 


where 
L t/2 


n,)? dn 


L 
liy+n, )(z+n,) dn 


F,*= 


ds fr (y+ny) dn 
o 
t/2 

— | as (z+n,) dn 


-t/2 


L 2 


| ds (yt ny) dir 


2? Cos ad 
Mott” COS 


—t/2 
L 
= I 2)| rue cos xls 
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(.\.59) 


(.\.60) 


The expressions given for the warping moments F,* and F,* may be transformed: 
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since the other two terms of the product give zero when integrated between the 
limits —t/2 and ¢,2. The definition of F, is given in Eq. A.20. Similarly, one 
may write 


Sin ads .  (.\.6rb) 


It should be remembered that in Eqs. A.60 and A.61 y and z are the co-ordinates 
of the median line measured from the centroid of the section. It may be seen 
that in the case of most thin-walled (t very small) sections the difference in the 
location of the centre of twist according to Eqs. .\.59 and according to Eqs. .\.22 
is insignificant. 

For the calculation of the shearing stresses associated with non-uniform warping 
one may assume that the component parallel to the median line is the same as 
that calculated in Sec. A.4 and given in Eq. A.24:— 


Ss 


except that the warping wy, corresponds to twisting about the centre of twist as 
defined by Eqs. A.59 rather than by Eqs. A.22. It should be noted that w, 
and not w,* figures in this equation. With this assumption the torque 7 of 
these shearing stresses with respect to the centre of twist is again given by 

L 

T=- | 

0 
On the other hand the shearing stresses according to Eq. A.24 have a non- 
vanishing resultant if the wall thickness of, the section is considered finite. This 
is seen when one calculates the components of the resultant in the same way 
as in Sec. A.4. One obtains 


s L 
L 
{ [ - | + 
0 


The first term in the twisted brackets vanishes for reasons stated in Sec. A.4. 
The second may be calculated from the first of Eqs. A.58 in conjunction with 
Eqs. A.55a and b. One has 


L t/2 
-t/2 
L L 1/2 
— [as | roan’ cos adn 
0 -t/2 
L L 


= i| — (t*/ 12) cos 


Consequently 
L 


Similarly one obtains 


Se= Eq!" (t/a tm sinads . (A\.62b) 
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The shearing stresses according to Eq. A.24 do not satisfy the equilibrium 
conditions of a small element cut out of the bar. The conditions may be satisfied, 
however, if one assumes that, simultaneously with the stresses according to 
Eq..\.24, shearing stresses 7, are acting in the direction of the perpendicular 
to the median line. From the theory of elasticity one has 

(Oo, + (07/08) + (07,,/dn) =0. 
This equation may be transformed with the aid of Eqs. A.55a, A.55b, and A.24 :— 
(We — lyn) — + (07, /dn) =o. 
Consequently, 
n 


fe) 
(n?/2) +f (x, 8). 
The unknown function f(a, s) can be determined from the requirement that the 
shearing stress 7, must vanish at the free surfaces n= —t/2 and n=t/2. Thus 
= (1/8) (4n? —t?) . (A.63) 
The y-component R, of the resultant of these shearing stresses is 
L t/2 L 
R, =f ds cos xdn= — (t*/12) cos zds .  (A.64a) 
o -t/2 
Similarly, the z-component R, of the resultant is 
(3/12) sin ads :  (A.64b) 
0 


It may be seen that the resultant according to Eq. A.64 just balances the 
resultant S according to Eqs. A.62. The moment T of the shear stresses 7, 
with respect to the centre of twist is 


L t/2 
[as 
-—t/2 


L 
2)) 
0 


This may be written in the form 


where 


Altogether the moment of all shearing stresses associated with non-uniform 
warping is 
where I’ is the quantity defined in Eq. A.28. Since I" is proportional to t, 
and I’, to ¢*, the latter may be neglected for most thin-walled sections. 
Exceptions to this rule are, for instance, the angle- and T-sections, for which 
I’ vanishes. 


Ag. THe WarpiInG oF Common AIrncRAFT SECTIONS. 


In the following the warping w and the warping constant I’ (and/or I’,) are 
given for a few sections common in aircraft construction. The sections are 


| 
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shown in Fig. A6. 


I- and T-sections, 


a defin 


disregarded, except in the cases of the angle- and T-sections, where I’ is zero 
for a rotation about the centre of twist. 

All the formule given are valid for a rotation about the centre of twist. 
the rotation is about some other axis, the value of w may be calculated with 


t 4 


{ 


ition cannot 


The co-ordinate s is measured along the median line of the 
section continuously from one free edge to the other, except in the case of the 
where such 
thickness ¢ is assumed to be small and the warping across the thickness is 


be maintained. 
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the aid of Eq. A.16 from the formule given below. Eq. A.16 may also be 
used, together with Eqs. A.18a and A.28, to develop a formula for T, for a 
rotation about any point P in the plane of the section :— 

(Yp— Yq) (2p —2q) Lyn t+ (Yp— Yq)? Len + (2p —2q)? - (.A.67) 
where |’, is the warping constant for a rotation about the centre of twist 7, 
yp and z, are the rectangular co-ordinates of the point P, y,; and z, those of 
the point T measured from the centroid C of the section, and the second moments 
of the ‘section [,,, [,,, and I,, are defined in Sec. A6. 

A. Channel Section.—The pertinent formule are given in Sec. A7. 

B. Z-Section.—The shear centre, or centre of twist, of the Z-section coincides 
with its centroid, and is in the middle of the side ) in Fig. A.6a. If the 
section is twisted about its centroid, the warping is 

w= —(h/2)b+C, when bSs S (b+h) (web) 
w= —(h/2)s+C, when oSs Sb (upper flange) (A.68) 
w=(h/2)(2b+h—s)+C, when (b+h) Ss S (2b+h) (lower flange) 

The warping constant is 

where 
k=(1/12)(b/h)* [4 +12 +9 (bh) + 2] /[8(b 
A.bgi 
+12(b/h)?+6(b/h)+1] 

C. Angle Section.—The shear centre, and consequently the centre of twist, is 
the intersection point of the median lines of sides a and b in Fig. A.ob. 
The warping is, therefore, zero everywhere, and so-is 

warping w,, however, does not vanish, and the warping constant 
1, = / 36) (a? + (.A.70) 
1). Flunged Angle Section.—The centre of twist T of the flanged angle section, 
shown in Fig. A.6c is on the axis of symmetry of the section at a distance 
e from the edge of the section :— 
e=(a/¥ 2) (b/a)? [3-2 (b/a)]/ [1 +3 (b/a)—3 (b/u)? +(b/a)?] (A.71) 
The warping of the section is given by the equations 
w=[at(e//2)] when oSs Sh 
w=(e//2)(s—a—b) when b Ss S 2a+b (A.72) 
w=[a+(e//2)] b—s)+(e/ 2) a when Ss S 


The warping constant is given by the equation 


where 
k=(1/6) (b/a)® [44+ 3 (b/a)]/[1+3 (b/a)—3 (b/ay + | (A.730) 
E. I-Section.—In the case of the I-section the co-ordinate s, is measured 
along the upper flange from s,=o to s,=h; the co-ordinate s, is 
measured along the lower flange from s,=o to s,=b. Both are shown 
S in Fig. A6d. The centre of twist coincides with the centroid. The 
: warping is given by the equations 


w=(1/4)h (b—2s,) upper flange 
w=o web ‘ (A.74) 

w=(1/4) h (b—2s8,) lower flange 

The warping constant is 
T=(1/24) (b/h)* th® : (.\.75) 
F. T-Section.—The case of the T-section is similar to that of the angle section. 
With the notation of Fig. A6e the warping constant is 

= (#9 / 36) (a* + (1/4) . (A.76) 
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Sus-ATomic Puysics. 
H. Dingle, Nelson's \eroscience Manual. T. Nelson and Sons, Ltd. 1942. 
5/- net. 

This book forms Volume II of a course of Physics for Aeronautical Students of 
which Volume I on Mechanical Physics has already been reviewed in the Journal. 

Like its predecessors in this excellent series, this is a first-class book. It covers 
in an introductory chapter the structure of the atom, and in the remainder of the 
book optics and electricity and magnetism. Professor Dingle not only writes well, 
but clothes the most difficult part of the book in common place analogy and 
language which go so far to make the student not only understand what the 
author is writing about, but also make him interested in what he is reading; and 
interest is ninety per cent. of the battle of educating. Nor does Professor Dingle 
shirk the absence or inadequacy of an explanation by a tangle of verbiage. 
Having spent over twenty pages describing in some detail the structure of the 
atom he adds: 

‘* The picture of the atom and its behaviour which we have given is the one 
which enables us to explain the largest number of facts without becoming too 
abstruse, The latest theories explain more but are not so easily grasped, present 
no imaginable picture of the atom, and are, moreover, liable to supersession by 
still more comprehensive theories as progress continues. The reader must under- 
stand that we do not say the atom is exactly as we have pictured it; it appears, 
in fact, that a precise picture is impossible to give. The value of the view we 
have adopted is that it gives a satisfactory and not misleading account of the facts 
which concern us in practical work. Where it fails we shall simply describe the 
facts without attempting to explain them.”’ 

Each chapter is followed by a series of useful exercises for the student. The 
book is excellently printed and illustrated and can be thoroughly recommended 
to all students taking the subject of physics in their examinations, whether aero- 
nautical or not. 


Atk NAVIGATION: British Empire Epirion (2nd Ed.). 
By P. V. H. Weems. McGraw Hill. 1942.  35/-. 

Any comment upon this work would be superfluous, since it has for many 
vears been established as a classic. The durable value of its contents is emphasised 
by the fact that practically nothing in the way of theory, methods or examples 
has been altered from the first edition. The changes brought about by five 
years of progress mainly concern improved types of navigation instruments. We 
have, for example, a description and photograph of the Astro Compass (pp. 
99-104); Drift Recorder, Mk. II; Course and Speed Calculator, Marks II and 
Ila; Navigational Computer, Mark III; and in a later chapter a description of 
the R.A.F. Bubble Sextant, with a much fuller account of methods of observation 
by Bubble Sextants and the efrors to which they are liable. Another useful 
addition is the method of finding fixed objects by Sun position lines (p. 366). 

To make room for this additional matter, very little has been omitted. Chapter 
XX, dealing with Navigational Equipment, has been entirely left out, presumably 
as being a mere catalogue of instruments that are described elsewhere ; while a 
further three pages have been saved by cutting out some notes on outstanding 
navigation flights, which, though interesting in themselves, are now ancient 
history. A minor, but significant alteration is the omission of the title ‘‘ Colonel ”’ 
before Lindbergh’s name—sic transit gloria mundi! 
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ASTRONOMICAL NAVIGATION. 
By W. M. Smart. Longmans, Green and Co. 1942. 5/-. 


This is the promised sequel to Professor Smart’s ‘* Introduction to Sea and 
Air Navigation ’’ and should be read in conjunction with it. The subject is 
handled in the lucid and practical manner which we have come to associate with 
the author. After a preliminary chapter giving the essential definitions and 
formule for spherical triangles, the book goes on to deal with the Celestial Sphere, 
the Astronomical Triangle, Greenwich Mean Time, Altitude Observations, the 
Astronomical Position Line, and the Air Astronomical Navigation Tables, while 
a concluding chapter is mainly concerned with Altitude Curves and other short- 
cut methods used by practising navigators. An Appendix reproduces some pages 
from the ‘* Air Almanac ’’ and the ‘‘ Astronomical Navigation Tables.’’ Each 
chapter contains a series of Exercises, to which answers are given at the end 
of the book. Star charts of the Northern and the Southern Heavens are printed, 
respectively, on the fly-leaves at the beginning and end. No space. has been 
wasted, but the quality of paper and printing is much better than in many 
wartime editions. In general, this is a book which should certainly be studied 
by all student navigators and might well be read as a ‘‘ refresher course ’’ by 
those who are more advanced. 


HANDBOOK. A TO AERONAUTICS, ATRCRAFT ENGINES AND SUPER- 
CHARGERS. (3rd Edition.) 
Flight Publishing Co., Ltd. 1942. 6/- net. 

The third edition of this most useful handbook has been considerably enlarged 
and a new section added on typical aircraft power units and supercharging. Many 
new illustrations, too, have been added, bringing the book completely up-to-date. 

There is not much doubt ‘* Flight Handbook ”’ will grow larger and more 
popular with each issue. For anyone with no knowledge of aeronautics, and 
eager to learn, this book will give him all he wants, in clear language, from the 
aerodynamics of the aeroplane to an explanation of the working principles of the 
instruments which are of such vital aid to the pilot. The new engine section is 
particularly well-illustrated with cut-away diagrams which mean so much to those 
who are becoming acquainted for the first time with engine construction. 


Gas TURBINES AND JET PROPULSION FOR AIRCRAFT. 
G. Geoffrey Smith, M.B.E. Flight Publishing Co., Ltd. 1943. 3/6. 

This is an important book. Mr. Geoffrey Smith has made a very thorough 
search for the work which has been done on gas turbines and jet propulsion for 
aircraft throughout the world. With great skill he has written his book in as 
simple terms as possible and has given to the young aircraft engineers, in whose 
hands the future of aircraft design lies, enough food for thought to keep them 
healthy for many years to come. 

‘* The idea of an aircraft being propelled without the aid of an airscrew strikes 
one at first as fantastic,’’ says the author. It might have the advantage of 
settling once and for all whether the term propeller or airscrew should be used. 
There is little doubt, however, that jet reaction has come to stay in some form 
or other and that it will play a major part in powering air transport of the future. 

The advantages of jet propulsion units or reciprocating engine units are such 
that intensive research and development are bound to be encouraged by the 
governments of the world to bring about the practical results which’ appear 
possible. 

With jet propulsion safety fuels can be used with direct application of the 
power produced. Construction of jet propulsion units is lighter and the power- 
weight ratio better than those of the reciprocating engine. The abolition of the 
airscrew will lead to lighter undercarriages, and units can be completely enclosed 
in the wings and fuselage, with an increase in aerodynamic efficiency. 
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There are disadvantages, but most of these are those of development, and no 
serious disadvantage from a theoretical point of view presents itself. Jet propul- 
sion, as the author rightly points out, is particularly suitable for aircraft operating 
at high-altitude levels, a fact of some importance for the progress of civil aviation. 

The coming of the jet will bring about new ideas in design, and aircraft of the 
future will make some of the best aircraft of to-day appear like ugly ducklings. 
A new era in power for the aeroplane is dawning. That does not mean to say 
the reciprocating engine is dead, or even dying, but it does mean that undoubtedly 
for the large, long-range, high-speed aircraft of the future jet propulsion offers 
more possibilities than the aircraft driven by airscrews. 

Mr. Geoffrey Smith has done well to collect together existing information and 
to point out the possible future of the jet. In this Journal was published the first 
article written by Whittle, whose name is so frequently mentioned in Mr. Geoffrey 
Smith’s pages, as long ago as 1931, on the turbo compressor. 

Germany and Italy have paid special attention to jet propulsion and its possi- 
bilities. ‘The only thing which puzzles the reviewer is why Italy gave so much 
publicity to the Caproni-Campini on its flight from Milan to Rome in November, 


1941. 


Air Pitot TRAINING. 


Bert A. Shields. McGraw Hill Book Co., London and New York. 1942. 
24/- net. 
This book is essentially intended to provide, as the author states in his preface, 
in a single volume all the information necessary to pass the written examinations 
for a private and a commercial pilot licence.’’ The licences are American licences 
and the whole book is written on the basis of American practice. 

There are four sections, the first on aircraft and theory of flight, the second on 
aircraft engines, the third on meteorology, and the last on air navigation. 

Mr. Shields is himself a licensed pilot with over 5,000 hours flying. He is also 
a licensed aeroplane and engine mechanic and an instructor in the subjects he 
has covered in his book. He particularly stresses that ‘*‘ great pains have been 
taken to avoid all unnecessary technical terms.” 

The author has, indeed, gone out of his way to explain in straightforward, 
simple language his subjects, buc it does not always follow that such language is 
the clearest and most accurate for examination purposes. The reviewer wonders 
what marks an examiner would give to the candidate who was asked ‘‘ What is 
gravity ? ’’ and receives the reply ‘* Gravity is the force that gave early pioneers 
their greatest difficulty.”’ 

But these things are a matter of teaching. There is much in this book which 
will prove useful and informative to those wishing to acquire a knowledge of the 
why and wherefore of flight. It has the merit, too, of being excellently printed 
and illustrated. 


ae 


THE AVENGERS. 
By Charles Graves. Hutchinson. 1942. 7/6. 
SEVEN Pibots. 
By Charles Graves. Hutchinson. 10943. 8/6. 
These two novels are, respectively, the second and third in the series initiated 
by ** The Thin Blue Line.’’ They belong to that curious type of fiction which 
used to be called ‘‘ roman a clef,’’ for which the recipe is simple: take a real 


man named (say) Donald Black, call him ‘*‘ Ronald White,’’ describe in detail 
some flying adventure of his, altering the time, place, and possibly the type 
of his machine, and insert somewhere in the book an official photograph taken of 
some similar incident. The result is no doubt highly diverting to pilots, and 
those in their immediate circle who can identify the originals of the characters 
and incidents described, and to all those youngsters whose ambition is to join the 
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R.A.F. The latter, when they have followed Mr. Graves’s heroes through three 
books, will certainly be equipped with an extensive vocabulary. Personally, the 
reviewer found some difficulty in distinguishing the seven (to him) quite unchar- 
acterised heroes, and could have put either book down at any point without any 
urge to find out what happened to them; but, since 30,0co copies of ‘* The Thin 
Blue Line *’ and 20,000 of ‘* The Avengers *’ have been sold, it appears that 


PLANES EXPLAINED. 
Roger Tennant. Argus Press. 1942. 1/- net. 


This is a book which is written to appeal to that ever-increasing public which 
wants to know why aeroplanes, apparently called upon to do the same job, have 
such different shapes and characteristics. Here one learns why wings are 
elliptical, tapered, thick or thin, and so on; why there are low-winged, mid- 
winged and high-winged monoplanes; why some aircraft have one rudder and 
others two; why some under-carriages are fixed and others retractable; and so 
on. The book is illustrated by just those examples which will be of such use to 
pilots, to members of the Royal Observer Corps and others who are very much 
concerned with identifying aircraft. There is not the slightest doubt if they know 
the why of a particular shape they will remember it far more easily than if they 
have a little more to go on than that it is so. 

This is a most useful, interesting and informative shillingsworth. 


ADVANCED FLYING. 
Captain Norman Macmillan, M.C., A.F.C. George Allen and Unwin, Ltd. 
1943. 2/6 net. 

Captain Norman Macmillan is not only a good pilot but he is also a good 
writer, two things which a spate of war books shows do not often go together. 
Captain Macmillan, indeed, has a long and growing list of books in his name, 
most of them worth while reading many times. His wide experience from fighting 
in the last war to test pilot and flying in many of the different parts of the world 
on many types of aircraft, enables him to write with some authority. 

In his first chapter he deals with the flying of fighters and the overwhelming 
importance of the care of his eyesight to the fighter pilot. Captain Macmillan 
points out the great advantages the modern fighter has over the fighter in the last 
war, but rightly insists that in the ultimate fight it is the pilot who sees more in 
less time who stands the best chance. 

The difference between flying multi-engined aircraft and single-engined are 
clearly explained in the second chapter, one which is naturally followed by a 
chapter on flying on one engine. Captain Macmillan’s description of the return 
of a Hampden from Wilhelmshaven on one engine is not only first-class, but will 
give many a newly-fledged pilot confidence when forced to fly under similar 
conditions. 

In Chapter IV the author explains what is meant by the safety factors of an 
aeroplane, pointing out how worth while it is to know the strength of the machine 
you are flying, ‘‘ and helps you to decide the problem of whether to bale out or 
to fly the crate home like ‘ a piece of lace.’ ’’ The Forces in Flight, the subject 
of Chapter V, naturally follows. Chapter VI is concerned with the effect of 
manceuvres both on the aeroplane and the pilot. 

Chapter VII covers Taking-off. As the author remarks, ‘‘ Every take-off 
varies.’’ Captain Macmillan explains how’certain games improvea pilot’s take- 


off abilities, and the reviewer heartily agrees with his comment. ‘* The static-ball 
game of golf is not a pilot's game.’’ Chapter VIII is concerned with the causes 
and dangers by the stall, and the last chapter with landings and under-carriages. 

This is a very good half-a-crown’s worth and even the most experienced pilot 
will get his half-a-crown back in some useful piece of knowledge. 
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An INTRODUCTION TO THE STUDY OF MAP PROJECTION. 
By James Mainwaring. Macmillan and Co. 1943. 5/-. 

It is several years since there has been published a good elementary text-book 
on map projection. The present work is therefore particularly welcome, at a 
time when so many young men are studying navigation in the A.T.C. or pre- 
entry classes for the R.A.F. It is divided into three parts, entitled respectively : 
The Problem of Map Projection, The Process of Map Projection, and The Projec- 
tions Applied. Of these, the second is, of course, the fullest and most important. 
Cylindrical, conical and zenithal projections are dealt with in turn, with a profusion 
of excellent diagrams and the strict minimum of essential mathematics. 
Particularly valuable are the comparative diagrams (Nos. 25, 26 and 27) summing 
up the scale differences between the various types of projection. Diagram 44 
and the hints on the page facing it will also be of great use to the student who 
may be called upon to recognise any given projection in an examination. The 
book includes a series of test questions which the student may work out for him- 
self from the text, or with the help of a teacher (answers are given only in a few 
numerical cases). 


CoastaL COMMAND AT War. 
By ‘* Squadron-Leader Tom Dudley Gordon.’’ Jarrolds, 1943. 7/06. 

By a coincidence this book appears simultaneously with the Air Ministry 
official publication entitled ‘‘ Coastal Command.’’ Many of the illustrations are 
the same—or obviously taken on the same occasion—and very beautiful they are. 
For sheer technical perfection it would be hard to beat the aerial views of the 
Norwegian fiords, the treeless hills of Iceland, or the floating pack-ice. The 
reader will be well advised first to run through the Air Ministry booklet, then 
to re-read the same incidents in the present book, where he will find the ‘Shuman 
interest ’’ angle which is necessarily left out of the official account. From 
internal evidence, the author is a professional writer: he has a dry and witty 
style, a nice sense of characterisation, and the ability to tell a good story. We 
particularly enjoyed his anecdotes of the Rear-Admiral Pilot Officer (page 31) ; 
of the ‘‘ mis-spelled pleasantry '’’ from the Hudson navigator (page 36), and of 
the German commander’s signal: ‘‘ What is ‘ scuttle’ please? ’”’ On the 
more serious side, we have vivid descriptions of the Altmark reconnaissance, the 
work of Coastal Command over Dunkirk, and its crowning exploit in hunting 
the Bismarck. 


ELEMENTARY ELECTRICITY FOR RADIO STUDENTS. 
By W. E. Flood. Longmans Green. 1942. 1/-. 

This addition to the ‘‘Service Manuals ”’ series will be welcomed by A.T.C. 
cadets and prospective wireless operators. Most text-books on electricity, even 
those describing themselves as ‘‘ elementary,’’ assume that the reader knows at 
least the ABC of the subject, and this frequently means that he is confronted with 


terms, quantities and circuit drawings which he finds puzzling. The present 
booklet assumes zero knowledge, and proceeds to impart the ABC in a pleasant 
and lucid way. Every symbol is identified, every term explained in detail, and 


every explanation is followed by a simple numerical example, worked out in full. 
Although the book is only 62 pages long, it contains an index, exercises, and a 
set of revision problems. 
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